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R.P. ADGAOITKA.R 
Ph.P. 

Department of Electrical Engineering 
Indian Institute of Pechnolog 7 , Kanpur 
April 1979 

DYNAMIC EQUIVALENTS POR POWER SYSTEM STUDIES 

Two major problems in the stability studies of preacnt 
day lar^ scale power systems are (i) dimensionality of the 
model and (ii) numerical stiffness of the system equations. 
The problem of dimensionality arises because of the large 
size of the interconnected power system and the recent trend 
t* represent the dynamics of a generating unit (synchronous 
machine and the associated controls) in greater detail. In 
addition to the dimensionality, the full scale dynamic model 
•f a large power system is associated with the problem of 
stiffness in the numerical integration because it accounts 
for both slow aad fast phenomena in the system. Hence the 
digital simulation of large systems become expensive in terms 
of computer time. Since the digital simulation for different 
loading levels, network configurations and disturbances is a 
routine feature in the power system planning and operation, 
the use of full scale models is uneconomical and undesirable. 
Therefore, simplified and reduced order models of power 
systems, which, apart from alleviating the difficulties 
mentioned above, describe the response of the system within 
engineering accuracy, are constructed. The problem of 



dynamio equivalence is concerned with deriving such reducedl 
order models. This thesis addresses itself to some aspects 
of this problem for transient (nonlinear model) and dynamic 
stability (linear model) studies. 

The concept of study system and external system 
(Undrill 1971), which is generally used in the development 
of an equivalent for transient stability studies, is 
familiar to power engineers. Thus the reduced model of a 
power system for transient stability studies consists of the 
full scale representation of the study system and an 
equivalent of the external system. Such a reduced order 
model (which we will refer as the equivalent of systan) is 
used to obtain the response of the study system for local 
disturbances as it is affected by the external system. Among 
the equivalents for transient stability studies, the principle 
of coherency based dynamic equivalence is likely to be widely 
used in power industry because of its simplicity, accuracy 
and the fact that physical correspondence with the components 
of a power system is retained in the equivalent. The basic 
steps in the development of a coherency based dynamic 
equivalent are (i) identification of the groups of coherent 
generators valid for a set of fault locations in a given study 
system and (ii) dynamic aggregation of each coherent group 
into an equivalent generating unit. 

Identification of coherent generators constitutes a 
key step in the construction of a coherency based dynamic 



equivalent. Since the knowledge of coherent groups valid for 
a set of fault locations in the given study system is a "basic 
requirement for constructing the equivalent, the problem is to 
identify the groups for multiple fault locations in the study 
system. The method for grouping the generators should be simple 
and efficient in terms of computer time. The reliability of 
results is of prime importance because there should not be any 
need to verify the results from the base case swing curves. 
Pormation of coherent groups on the basis of base case transient 
stability studies results in prohibitively large computer times. 
Many of the existing methods of coherency analysis (lee, 
Spadling, Voropai and others) do not reflect the maximum 
angular excursion between the generators during the transient 
period and hence are empirical. The method of linear simulation 
proposed by Podmore uses simplified swing curves and has been 
proved to be reliable by extensive testing. However, the 
method requires storage and comparison of swing curves, 
Furthermore the linear simulation has to be repeated for 
every change in the fault location. Therefore, there is a 
need for a direct criterion for coherency analysis which is 
suitable for multiple fault locations and does not resort to 
storage and comparison of swing curves. 

Over the past few years, although the equivalising 
techniques have been perfected, the basis for determining the 
levels of complexity in modelling the genererbtars (before an 
equivalent is developed) has remained largely heuristic. For 



example, it has been a common practice to represent the 
generators avay from the fault location with lesser detail 
than those in the vicinity of it, thereby ignoring completely 
the electromechanical effects during the transient period. In 
the li erature, the admittance and reflection distance 
measures (Lee and Schweppe) are proposed for the above 
purpose. While the admittance distance is a s tat i c measure, 
the reflection distance partially reflects the electro- 
mechanical effects during the transient period. Therefore 
there is a need for an improved criterion which demarcates 
a power system into different regions for the purpose of 
modelling the generators. 

The existing techniques for simplifying the linear 
dynamic models of power systems for dynamic stability studies 
are mainly based on modal analysis (Kuppu2?a3ulu, Altalib, 

Van ITess and others) i.e. retaining the dominant modes and 
eliminating insignificantly excited modes. This, however, 
involves the computation of eigenvectors and reciprocal 
basis vectors, which becomes unwieldy if the size of the 
system matrix is large. Therefore, the time required to 
compute the reduced model may prove to be prohibitive unless 
the equivalent is used for multiple stability studies. The 
state variable grouping technique (Kuppurajulu) does not 
require these computations but the method requires apriori 
knowledge of the behaviour of the state variables. Thus 
there exists a scope for developing simpler techniques with 



a rigorous mathematical basis for simplifying the linear 
models of power systems for dynamic stability studies. 

This thesis presents solution to some of the above 
mentioned problems. The scope of thesis is broadly divided 
into tw. parts, 

(i) Identification of coherent generators and decomposition 
of power systems using a new approach, and construction 
of coherency based dynamic equivalents. 

(ii) Application of the singular perturbation theory to 
obtain simplified models and simulate the power system 
for stability studies. 

A chapterwise summary of the work done in this thesis 
is given below. 

The first chapter includes introduction to the problem, 
a brief review of the existing methods and the scope of the 
thesis. 

The second chapter deals with a new method for 
identifying the coherent generators in power systems. The 
method j.s analytical and uses the concept of coherency index 
for grouping the coherent generators. The linearised swing 
equations of the system are cast into the usual state space 
form after selecting the generator nearest to the fault 
location as reference generator. Effect of fault is taken 
into account in terns of an input vector for the linearized 
system as done by Podraore. A closed form solution of the 
system equations is obtained in terms of the decoupled modes. 



The extent to which each mode is excited is determined hy the 
fault condition under consideration. A new concept of coherency 
index, which is a measure of the maximum angular excursion 
"between two machines during the transient period, is introduced. 
The col: srencj?- index exploits the nature of the closed form 
solution in terms of the eigenvectors, reciprocal "basis vectors 
and the state of the system at the time of fault clearing. 
Therefore the coherency index between two machines serves as a 
logical criterion to determine the coherency between the 
machines. "When it is required to determine the groups for 
different fault location in the study system, the eigenvectors 
and reciprocal basis vectors corresponding to a new reference 
generator are not recomputed. Instead they are derived from 
the existing ones by using a similarity transformation. The 
use of similarity transfomation eliminates the repetitive 
computations for different fault locations. In a similar 
manner the method can be easily adapted for the change in the 
geographical configuration of the study system itself. Two 
power system examples are presented to validate the theory in 
practical systems. 

In the third chapter, the problems of decomposition of 
power systems and construction of coherency based dynamic 
equivalents are taken up. The decomposition of a power system 
into different regions is primarily aimed at determining the 
relative degree of complexity in modelling the generators 
before an equivalent is developed. A new concept of 



electromechanical distance, which is used as the basis for 
decomposition of a power system, is presented. The concept 
is similar to coherency index except that it is measured 
between a generator in the external system and the reference 
generator in the study system. The decomposition of the power 
system is linked with coherency analysis and the possibility 
of grouping the generators in each region separately is 
established. The equivalising methods for combining the 
terminal buses of coherent generators and the aggregation of 
coherent units are described. A case study of a typical power 
system in India illustrating the principle of coherency based 
equivalents is presented. 

• The rest of the thesis is devoted to the 
application of singular perturbation theory to the problems 
of dynamic simplification and simulation in power system 
stability studies. The primary motivation for using the 
singular perturbation technique is that, in this technique, 
the slow and the fast subsystems are solved separately and in 
two different time scales. The technique first reduces the 
model by neglecting the fast phenomena in the system. It then 

improves the approximation by reintroducing their effects as i 

i 

’boundary layer correction’ calculated in different time scale. 
IPurther improvements in the solution are possible .by evaluatang i 
the subsequent terms in the asymptotic expansion. These 

j 

features of the singular perturbation theory make it very 
suitable for application in large scale systems. The technique 
has been extensively applied in optimal control and linear 



regulator problems by Kokotovtc and others. It has also been 
used in other branches of engineering. However, so far, the 
technique has not been used in the stability studies of large 
scale power systems for model simplification anfl subsequent 
Simula "t ton. 

The fourth chapter describes the application of singular 
perturbation theory to power system dynamic stability studies. 
Pirst a state space model of a multimachine power system for 
dynamic stability studies using the detailed representation of 
the dynamics of the generating units is obtained. These system 
equations are transformed into the singular perturbation form 
by properly identifying the perturbation parameter 'e ’ in small 
rotating masses, small time constants and large gains etc. The 
use of 'e' in practical systems is symbolic and represents the 
presence of fast as well as slow subsystems. The presence of 
these subsystems in a given model is identified by a criterion, , 
which examines the eigenvalues of the system matrix. The i 

degenerate system i.e. setting e = 0 in the system equations 
gives +he simplified model of the original system. In the 
simplified model, the identity of the corresponding state : 

variables of the original system is retained. The response ; 

of the degenerate system is further improved by adding the 
terms corresponding to 'e’ in the asymptotic expansion to 
obtain a first order solution. ^ single machine infinite bus i 
and a three machine example is presented to illustrate the 
theory. It is also established that the first order solution 



is a good approximation to the response oi the original system 
and can be obtained with far less computer time as compared to 
simulating the original system. 

Chapter five deals with the singular perturbation 
approach to the transient stability studies in power systems. 
The technique uses the method of asymptotic expansion as 
applied to nonlinear differential equations (initial value 
problem). The relevant theory is described and illustrated 
with the help of a numerical example. The difficulties 
encountered with nonlinear models of power systems are 
clearly brought out and explained. 

A summary of the contributions in the thesis and the 
scope for further research in this field are given in the 
concluding chapter. 



OHAPTER 1 


IFTROnJCTIOlT 


1 . 1 STABILITY STUDIES 

Stability considerations have been recognised as an 
essential part of pc^er system planning and operation for a 
long time. In the planning stage, routine stability calcula- 
tions are made to ensure the stability of the system with 
proposed transmission and generator configurations. The 
stability studies are also necessary for coordination of 
system protection and design of control systems. In the 
’on line’ operation of power systems, quick stability 
calculations at regular intervals enable the operator in an 
energy control center to take corrective actions so as to 
ensure security of the system against likely contingencies. 

A power system engineer is primarily concerned with two 
types of stability studies namely (i) transient stability 
and (ii) dynamic stability. While the former is associated 
•with large and aperiodic disturbances and uses nonlinear 
mathematical model, the latter is concerned with small signal 
analysis around an operating point. The disturbances consi- 
dered for dynamic stability studies are continuous and small 
in magnitudes so that the linearised mathematical model is 
adequate. 

1 . 2 STATEMENT OP THE PROBLEM 

The mathematical model of a power system for stability 
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representing the dynamics of the generating units (synchronous 
machines and their controls) and a set of algebraic equations 
describing the interconnections through the transmission network. 
Ihe major problems encountered in the dynamic simulation are 
(i) dimensionality of the model and (ii) 'stiffness' of the 
system equations. The problem of dimensionality arises due to 
the large size of a power system and is further aggravated by 
increased interconnections and the recent trend to represent 
the dynamics of the generating units in greater detail. It is 
highly uneconomical and undesirable to simulate a large sized 
model on a digital computer for different loading conditions, 
system configurations and disturbances. The other problem of 
numerical 'stiffness' of system equations arises because the 
system matrix has eigenvalues with very small as well as very 
large magnitudes. This necessitates use of prohibitively 
small time step size in numerical integration and hence the 
dynamic simulation becomes expensive in terms of the computer 
time. In order to reduce the computations and storage require- 
ments without sacrifice in engineering accuracy, one resorts to 
simplified and reduced order models. The problem of dynamic 
equivalence is concerned with deriving such reduced order 
models. Some of the important aspects of the problem are 
(i) time of computation required for constructing the 
equivalent, (ii) an appropriate compromise between reduction 
in size and loss of accuracy and (iii) physical correspondence 
of the reduced system with the components of the original 
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system. The last aspect makes it possible to use the conven- 
tioual stability programs, which have beoi perfected over the 
years, in simulating the reduced system without any modifi- 
cation. 

This thesis addresses itself to this problem of 
dynamic equivalence and simulation of power systems for 
transient and dynamic stability studies through reduced order 
models. The first part of the thesis is concerned with 
coherency based dynamic equivalents for transient stability 
studies. The following problems from this area are 
investigated, 

(i) Identification of coherent generators in large scale 
power systems based on eigenvectors and reciprocal 
basis vectors of the linearized system. 

(ii) Decomposition of power systems into regions using 
techniques developed for (i) to determine the 
relative degree of complexity in modelling the 
generators. 

(iii) Construction of the dynamic equivalent. 

The secrond half of the thesis deals with the 
application of singular perturbation theory to the problem 
of dynamic equivalence and simulation of power systems for 
transient and dynamic stability studies. The problems 
investigated are the following. 

(i) Modelling of power systems in the singular perturbation 
form for dynamic and transient stability studies. 
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(ii) Development of the reduced order model and the 

subsequent improvement in the response using singular 
perturbation theory. 

We now briefly review the pertinent literature in these 
areas before outlining the objectives and scope of the thesis. 

1 . 3 THE STATE OE ART 

1.3.1 Dynamic Equivalents for Transient Stability Studies: 

In the development of dynamic equivalents for transient 
stability studies, a power system is generally divided into 
two parts namely (i) the study system and (ii) the external 
system. The study system is that part of the system in which 
the system behaviour is of direct interest and where a set of 
disturbances are considered. The part excluding the study 
system is the external system, which is equivalised. Thus the 
reduced model of a power system for transient stability studies 
consists of the full scale representation of the study system 
and an equivalent of the external system. Such a model (which 
we will refer to as the equivalent of the system) gives the 
response of the study system for the disturbances in the study 
system. If the external system has been properly equivalised, 
this response should agree within engineering accuracy with 
that of the original full scale system. 

The existing techniques for constructing the 
equivalents for transient stability studae^ are bro.adly 
classified ini:© three groups. 
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(ii) Equivalents based on modal analysis [4-6], 

(iii) Equivalents based on coherency [?-14]. 

In (i), the generation and loads in the external system 
are equivalised based on the concept of distribution factors. 
However the resulting reduced order model may not preserve the 
dynamic characteristics of the original system to sufficient 
accuracy because the generators are equivalised irrespective 
of their different dynamic responses. The method of alloca- 
ting the inertias to the equivalent generators l3l is heuristic. 
Therefore these equivalents give very approximate results. In 
a recent paper [ 72 1, new distribution factors are defined based 
on the concept of moment. Here the dynamic effects of inertias, 
generation and loads in the external system are transferred to 
the boundary buses using the laws of mechanics. Por this purpose, 
the machine inertias, loads and generation at the respective 
buses are considered as forces and the transfer impedance 
between the two buses is considered as the distance. 

In the modal equivalent method [4], the external system, 
where the disturbance propagated is assumed to be small, is 
modelled by linearised differential equations and the study 
system by the usual nonlinear differential equations. The 
external system matrix is reduced by retaining only the 
dominant modes. A mode is neglected [ 5 ] if (i) it decays 
fast as compared to other modes or (ii) it is insignificantly 
excited. Modal equivalents need computation of eigenvectors 
and reciprocal basis vectors, which becomes unwieldy if the 
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size of system matrix is large. A recent comprehensive testing 
of the modal equivalent method 161 indicates the suitability of 
the technique only v/hen a number of faults in the study system 
are analysed using the same equivalent. The reduction of 
external system matrix needs a better understanding of mode 
selection process and any error in the judgement results in 
degradation of the performance of the modal equivalent. Further- 
more, the reduced model does not retain the physical corres- 
pondence with the components of a power system. 

In the coherency based equivalents 17-141, a group of 

generators having similar dynamic responses for a fault in 

the study system is replaced by an equivalent generator. 

References [7-9 1 describe algorithms for obtaining a reduced 

set of transient stability equations based on the concept of 

coherency. However, the reduced equations derived in these 

references do not reflect physical significance in terms of 

the normal power system components. A significant contribution 

in the area of coherency based dynamic equivalent was made by 

the research group of Systems Control Inc. [10-131. Their 

approach of constructing the equivalent is practical and is 

likely to be widely used in power industry. The basic steps 
¥ 

consist of the following. 

(i) Definition of the study system. 

(ii) Identification of groups of coherent generators in 
the external system whi'-^h are valid for a set of 
disturbances in the study system. 



(iii) Reduction of coherent generator buses. 

(iv) Reduction of load buses. 

(v) Dynamic aggregation of coherent generating units. 
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The algorithms of identifying the groups of coherent 
generators reported in [12] are covered in Section 1.3.2. In 
the reduction of coherent generator buses, the terminal buses 
of the generators in a coherent group are replaced by an 
eijuivalent bus rather than combining the internal buses of the 
generators [7,8 1. Such an approach decouples the problem of 
network reduction from that of dynamic ag=yregation of the 
generating units. The method used for reduction of load buses 
is selected based on linear or nonlinear representation of 
loads. A coherent group of generators is represented by an 
equivalent terminal bus and a generating unit connected to it. 
The procedure of determining the parameters of the equivalent 
generating unit from the parameters of the individual units 
in a coherent group is termed as dynamic aggregation of the 
generating units [12,151.. The components of generating units 
in a group such as synchronous machines, ejccitation-voltage 
regulator systems, turbine-governor systems etc. are aggregated 
separately. The technique for determining the parameters of 
an equivalent component is to numerically adjust the para- 
meters so as to obtain minimal error between the frequency 
response of the transfer function of the equivalent component 
and the sum of frequency responses of the individual transfer 
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functions. The error to be minimised is the difference of 
squares in the two frequency responses at specified discrete 
frequencies . 

The approach of Systems Oontrol Inc. [12] has the 
following features. 

(i) The development of equivalent is systematic and 
general, accommodating different types of synchronous machines, 
their control equipments and the nonlinear characteristics of 
loads. 

(ii) The equivalent is composed of normal power system 
components and, therefore, is compatible with the existing 
transient stability programs. 

(iii) The different tests on large scale power systems [12] 
demonstrate the capability of the equivalent to handle the 

transient stability simulaiii-OTi of large systems within a reason- 
able time. 

Other references on the equivalents for transient 
stability studies are given in [14-21]. 

1.3.2 Coherency Analysis: 

Identification of coherent generators in the external 
system valid for a disturbance/a set of disturbances in the 
study system is a key step in the development of coherency 
based dynamic equivalents. In the earlier stages of its 
development, the base case transient stability swing curves 
were compared to classify the coherent generators. However, 
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the numerical integration of nonlinear equations is time 
consuming and hence undesirable for the purpose of coherency 
analysis. Hence a need arose for more efficient and simple 
algorithms for coherency analysis and many researchers have 
contributed to this area. 


(a) Pattern Recognition Method : Lee and Sohweppe [22], based 

on extensive analog simulations, presented three features 

which are central to the phenomenon of coherency. The features 

are as follows: 

Max(a.) - Min (a.) 
i ^ i ^ 


^1 = 


Ma!x (a^j^) 


( 1 . 1 ) 


where a^i^ is the acceleration of generator 'i’ at t=o' 


r2 = 


Min r Min (y . • )1 
ieO Lie0..i^i ^ ^ f 


Max 

ieC 


r Max (y., ) 


( 1 . 2 ) 


where C = set of machines considered for coherency, 

4> = set of inner circle (study system) generators, 
y. . = transfer admittance between the generators 'i’ and 


’ 0 ' 


r, = 


Max(M. ) - Min(M. ) 
1 ^ i ^ 

Max(M. ) 
i ^ 


(1.3) 


The coefficient compares the degree of difference 
between the accelerations at t=o''‘ of the machines considered 
for coherency. r 2 measures the electrical coupling between 
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the machines considered for coherency relative to the machines 
in the study system. The feature r^ defines the similarity 
of inertias. A measure of the degree of coherency is defined 
by the coherency factor as 


y = Max 
t 

The coherency factor defines the allowable deviation in the 
relative rotor angle among coherent machines as a percentage 
of the maximum deviation in the transient period. 


Max 


A 6^(t) - A6^(t) 
Max( ^ 


o < t < 1 sec. 


(1.4) 


The coherency criteria for a 10^ coherency factor 'y' 
is defined as 


r.| ^ 30 /!^ , ^2 — ^ 


(1.5) 


As it is clear from the definitions, the features r^ , 
r 2 and r^ are defined by the system parameters and the 
acceleration ^wer at t=o''’. These features do not reflect the 
maximum deviation between the machines during the transient 
period. Therefore the method is empirical and lacks the 
required consistency in results [12], 

(b) Technique using Singular Point : Spadling et al.[ 23J 
presented a technique for identifying coherent generators 
by making use of the singular point concept. In this 
technique, the rotor angles of the generators are compared 
at both the initial stable operating point and the nearest 
post-fault unstable equilibrium point corresponding to the 
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expected mode of instability for a given fault. Pig. 1 .1 
shows the region of attraction around the origin, which is 
the post-fault stable operating point. The pre-fault stable 
operating point (initial conditions) is denoted by 6 ^ and the 

s 

nearest unstable equilibrium (singular) point is denoted by £ . 
The authors claim that if the fault is cleared at 12? !•©• 
critical clearing time, the resulting trajectory is maximal 
and approaches the nearest unstable equilibrium point. However 
this fact is not observed in the base case transient stability 
studies because the singular point under consideration is 
unstable. Two generators 'i’ and ’j' are classified as 
coherent if 

|6?. - 6?.l < e (1.6) 

13 13 - 

where e is the tolerance. 

Th.ei computaihL'onal. effort required by this method is 
very small i.e. the solution of a set of(n-l) nonlinear 
algebraic equations for a n machine system to obtain the 
required singular point i®. For this purpose, the mode of 
instability for a given fault has to be assumed apriori. This, 
in general, is a very difficult problem and depends on type, 
location of fault and system loading conditions. The 
QPi 02 ^ j_Qn (1 .6) verifies the features of coherency only at 
two points of the trajectory. Generators which are geographi- 
cally remote may satisfy the criterion (1.6). Therefore the 
criterion is to be backed up by the concept of admittance 
distance [221. 



12 





^2 = Critical fault clearing time 


Fig.M Trajectories in a region of attraction 
around the post -fault stable 
operating point 
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(c) L inear Simulation Method ; The method of linear simulation 
[ 243 , proposed by Podmore, essentially determines the groups of 
coherent generators by comparing the swing curves obtained from 
the simplified linear model of the system. An important aspect 
of this method is the manner in which a fault condition is 
reflected. The response of a faulted system is approximated by 
the response of the unfaulted system when each generator is 
subjected to a step input. The magnitude of the step input 
to a generator is equal to the accelerating power of the 
generator at t=o'*’ and duration of the input is equal to the 
clearing time of fault. 

Considering the classical models for the generators, 
the linearised swing equations for mechanical motion of the 
ith synchronous generator in a ’n' machine system is given by 


d ( Ao) . ) 

M. = AP . - AP . - L. Ao), 

1 dt mi gi 1 1 


the d.c. load flow equations are 


(1.7) 

i=1 ,2, . . . ,n (1.3) 

(1.9) 


where P_ and d = the powers and rotor angles at generator 

~6 ~ 

internal buses, 

P^ and _0 = the powers and voltages angles at load buses, 





Jl 


- 





[ 

A^ 

i 


The elements of the Jacobian H are derived from the appropriate 
partial derivatives or, if the flat voltage profile is assumed, 
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it may 'be derived directly from the [B] matrix of the fast 
decoupled load flow [251. The different types of disturbances 
are reflected in depending upon the electrical 

fault, loss of load or loss of generation. The swing curves 
are obtained by integrating the equations (1.7) and (1.3) by 
trapezoidal rule [261 and solving the network equations (1.9) 
at the end of each integration step. The groups of coherent 
generators are determined by processing the swing curves in a 
clustering algorithm t12l. 

The modelling of disturbances proposed in this method is 
valid for typical fault clearing times and the swing curves 
obtained from the linear model are shown to preserve the 
coherency property in the system. Therefore, the method is 
simple and efficient. The extensive testing of the method on 
different practical power systems [121 has proved the reliabi- 
lity of the results. However, the simulation has to be 
repeated for every change in the fault location and the 
method needs storage and comparison of swing curves. 

(d) Other Methods : The use of Laplace transform and Fourier 
transform have also been attempted in [12, 301. However these 
approaches have not been found computationally attractive as 
compared to linear simulation. Schlueter et al.[27] derived 
an expression for r.m.s. coherency measure [ 11 ] using the 
linearised model [24 1. Their method uses probabilistic model 
of power system disturbances and the concept of modes in 
deriving the expression for r.m.s. coherency measure. 
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Pai and Farayana [141 viewed the coherency between the 
generators from the energy concept. In their approach, two 
generators 'i' and ' j' are classified as coherent if the 
potential energy component associated with this pair is 
zero during the transient period. Mathematically the 
criterion is expressed as 




Voropai [231 derives certain coherency coefficients 
from the following ejuation for a pair of machines. 


h: 


P. . - P. . Sin <S. . - A. . 
ID ID ID ID 


( 1 . 11 ) 


where P. P. ^ and A. . are assumed as constant. The concept 
of equal area criterion is used to derive these coefficients. 
An examination of the above equation shows that the inter- 
connections of the machines ’i' and ' d' v/ith the rest of the 
system are completely ignored. Therefore, such decoupled 
equations are not valid in multimachine systems. The ’ct* 
and '6' coefficients for coherency analysis derived by 
Varwandkar [291 are also constrained by a similar drawback. 


The above survey of coherency analysis indicates that 
there are many methods of coherency analysis in the literature. 
However, many of them lack the mathematical basis, consistency 
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of results or they are uot tested for large scale systems. The 
most promising technique, however, is that of linear simulation 
{ 12 , 24 ]. The approximations involved in this method have been 
validated on large scale power systems. The method is simple 
and efficient. However, it requires the computation and 
comparison of swing curves and the process has to be repeated 
as the fault is moved in the study system. The same is true if 
the geographical configuration of the study system is changed. 

The method presented in this thesis overcomes some of these 
drawbacks. Our approach using the linearised model, develops 
coherency indices derived from the eigenvectors and reciprocal 
basis vectors for determining the groups of coherent generators. 
The method is shown to be useful when different fault locations 
are considered or the study system gets changed. Its demonstra- 
tion for a large realistic systems establishes its utility for 
practical applications. 

1.3.3 Simplified Models for Hynamic Stability Studies: 

The techniques for simplifying the linear models of 
power systems used in dynamic stability studies are mainly 
based on modal methods [31-331. In the modal methods, the 
dominant modes are retained in the reduced model and the 
insignificantly excited modes are neglected. Kuppurajulu and 
E]langovan[39 1 derived three different reduced models for initial, 
intermediate and final responses using the eigenvalue grouping 
technique. Davison’s method [31] has been applied to simplify 
the linear dynamic models of power systems by Altalib and 



17 


Xlrause £411. They have used a decentralised approach in the 
sense that, first the reduced order model of each subsystem 
such as exciter-voltage regulator system, turbine governor 
system etc. is reduced separately and these individual reduced 
order models are combined to obtain the reduced model of the 
system. Radhakrishna et al.[42] proposed the decomposition of 
full scale system model into different subsystems based on the 
criterion for weakly coupled systems given in [43J. Each of 
these subsystems is reduced individually and the reduced sub- 
systems are combined to obtain the equivalent of the original 
system. These decomposition techniques, however, need to be 
tested for large power systems to establish their practical 
utility . 

In the topological method presented by Van "Mess [ 441, 
the block diagram and the interconnections that are to exist 
in the reduced order model are apriori specified. Such a 
specification of the topology determines the nonzero elements 
oflAlmatrix. As in the modal methods, the eigenvalues of the 
reduced system are chosen from the eigenvalues of the original 
system. Some elements of eigenvectors are also chosen but 
there must be some relaxation of eigenvectors in order that 
the required topology can be realised. Then the values of 
nonzero elements of[ A] matrix, which give the required eigen- 
values and eigenvectors , are computed. The method is not simple 
to use in practice because in some cases the solution may not 
exist and it is difficult to determine the parameters 
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that will give the desired eigenvalues and eigenvectors. In 
the suhsequent paper [451 the author used the optimisation 
techniques to determine the parameters of the specified 
reduced model. The difference between the specified and the 
actual eigenvectors is minimised subject to the constraints 
that the re dured matrix has the specified eigenvalues. 


The methods for simplifying linear models of power 
systems, which do not resort to the eignveotors and reciprocal 
basis vectors, are the classical method- [441 and state variabl 
grouping technique [391. In the classical method, the dynamic 
elements such as time constants, rotating masses etc., whose 
effect on the system response is considered to be negligible, 
are deleted. The classical method is simple but gives the 
least a ccurate representation of the original system. In the 
state variable grouping technique, the state variables X of 
the original system 

X =[Al X + IBJ u ; X(o) = X® (1.12) 


are grouped as 
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(1.15) 


The groups X^ and X 2 associated with large and small 

time constants in the system respectively. The state 
variables X^ (which vary slowly as compared to X 2 ) are assumed 
to be constant during the initial period. Then the transient 
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response in the initial period is represented hy the simplified 
model 

Xi = 

X 2 =[A22lX2 +^^2^11 ( 1 . 14 ) 

During the final stages of the transient response, the fast 
traasients would have decayed and X 2 = £• Hence the system 
model for the final stages of the transient response is given 

X^ =tA*]X^ +[B*lu (1.15) 

-1 -1 
where [A*] = [A^ ^ -A^ 2 ^22 "^21^ [B*] = [B^-A^2 -^22 ®2^ 

The method is simple and intuitively correct. But the 
grouping of state variables is not backed up by mathematical 
criterion and one has to have apriori knowledge of the behaviour 
of the state variables. The boundary separating the periods of 
initial and final responses is not clearly defined. Later we 
will show that the models (1.15) and (1.14) can be derived from 
the singular perturbation approach. 

An inspection of the preceding methods for the simpli- 
fication of linear models of power systems brings out the 
following points: 

(i) Modal techniques, though mathematically sound, need a 
clear understanding of mode selection process and computation 
of eigenvectors and reciprocal basis vectors. 
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(ii) State variable grouping technique is efficient but 
needs the apriori knowledge of the behaviour of the state 
variables. Therefore the validity of results largely depends 
upon the human judgement. 

The above review underlines the need for a distinct 
and simple approach to the model simplification and simulation 
problem in dynamic stability studies. In the approach presen- 
ted in this thesis, the response of the full scale model is 
considered as a perturbation of the response of the reduced 
system. Singular perturbation theory is well known in the 
mathematical literature. In the next section, we briefly 
review the important aspects and the applications of the 
singular perturbation theory pertinent to the power system 
problems. 

1.5.4 On the Singular Perturbation Technique to Power System 
Problem: 

Two important features of the singular perturbation 
technique which make it suitable for large scale power systems 
are the following: 

(i) firstly the technique reduces the size of the given 
model by neglecting the fast phenomena in the system. It then 
reintroduces the effect of the fast phenomena as a correction 
term known as the boundary layer correction. 

(ii) The fast and slow subsystems are solved separately 
and in different time scales, thereby alleviating the stiff- 
ness problem associated with the original system of equations. 
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An excellent review of the singular perturbations and 
order reduction in control theory is given in Ref. [46], The 
paper also gives a thorough list of references on the theory 
and applications of singular perturbation technique. A few 
important references which o.re useful for the application of 
the technique in power system studies, are given in C 46-603 . 
Although the technique has been extensively used for model 
simplification in control theory and other branches of 
engineering, so far it has not been used in the studies of 
large scale power systems. A recent paper by Ihow et al.[57] 
deals with the singular perturbation analysis of systems with 
sustained high frequency oscillations. In this paper the 
theory is illustrated by a three machine power system example, 
where the high frequency oscillations arise between two 
coherent machines. In the technique for model simplification 
of power systems presented in this thesis, we have basically 
used the solution of the initial value problem in terms of 
asymptotic expansions given by O'Malley r503. 

1.4 SCOPE AND OUTLINE OP THE THESIS 

With the preceding background of the work done in the 
area of dynamic equivalents and model simplification for power 
system stability studies, we now present the main objectives 
of the thesis and a chapterwise summary. 

The main objectives of the thesis are as follows: 

(i) Development of an effective criterion for coherency 
recognition which does not resort to simulation and comparison 
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of swing curves and is suitable for determining the coherent 
groups valid for a set of fault locations in the study system. 

(ii) Decomposition of the power system hy means of a 
criterion into different regions for determining the relative 
degree of details in modelling of the generators. This criterion 
is consistent with the criterion for coherency analysis. 

(iii) Development of a technique based on the singular per- 
, turbation theory for the model simplification and simulation 

in transient and dynamic stability studies. 

A chapterwise summary of the thesis is given below. 

Chapter 2 is concerned with the problem of identifica- 
tion of coherent generators. The problem is introduced by 
enumerating the desirable features of a method for coherency 
analysis. The assumptions of linearity, classical models for 
generators and a typical modelling of electrical disturbances 
[24] are made in the formulation of the mathematical model for 
coherency analysis. Selecting a generator near the fault 
location as reference, the linearised swing equations are 
cast into the usual state space form as 

X = [A] X + [Bl u (1 .17) 

The forcing function u, which is the vector of accelerating 
powers of the generators at t=o''', reflects the fault under 
consideration. The vector X and the matrices [A] and [B] are 
defined separately for the cases of uniform and non-uniform 
damping in the system. Assuming that the matrix [A] has all 
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distinct eigenvalues, the closed form solution of (1.17) is 
obtained in terms of the decoupled modes. A new concept of 
coherency index, which is a bound on the maximum deviation in 
the rotor angle between two machines during the transient 
period, is introduced. The expression for coherency index is 
derived from the closed form solution of (1.17). The coherency 
index 'h- is used as a criterion for identifying the 

X J 

coherency between the generators 'i' and 'j'. A systematic 
procedure for grouping the coherent generators is presented. 

In order to determine the coherent groups for a different 
fault location in the same study system or in a different 
study system, a change of the reference generator may be neces- 
sary. In such cases, the eigenvectors and the reciprocal basis 
vectors corresponding to a new reference generator are not 
recomputed. Instead they are derived from the existing ones 
by using a similarity transformation matrix. The similarity 
transformation matrix can be easily obtained by inspecting 
the state variables with respect to the new reference 
generator. This feature of the method eliminates repetitive 
computations of eigenvectors and reciprocal basis vectors, 
when different fault locations are considered. The technique 
of grouping the generators and the use of similarity trans- 
formation is validated on typical power system examples. 

The third chapter deals with the decomposition of the 
power system and linking of the concept with the grouping of 
coherent generators. The concept of electromechanical 



24 


distance which is also a hound on the maximum deviation of a 
generator with respect to the reference generator in the study 
system, is used as a criterion for decomposing the external 
system into three or more regions , The study system, hy 
definition, constitutes the innermost region. There are two 
motivations behind the decomposition of the external system. 

(i) It helps in determining the relative degree of detail 
in modelling of the generators. 

(ii) The concept of electromechanical distance is used to 
delineate the boundaries between the regions in such a manner 
that a generator in one region is not coherent with a generator 
in another region. This simplifies the procedure for formation 
of coherent groups. 

The theory of such a unified approach for decomposition 
and coherency analysis is presented and validated on a power 
system example. Chapter 3 also contains a case study for a 
typical power system illustrating the principle of coherency 
based dynamic equivalents. 

'Chapter 4 addresses itself to the problem of model 
simplification and simulation of the power system for dynamic 
stability studies using singular perturbation technique, 
firstly the singular perturbation theory as applied to the 
initial value problem [501 is explained. The important 
features of the technique are brought out and the motivation 
behind the application of the technique to the problem of 
dynamic stability studies is stated. A systematic procedure 
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for formulating a state space model of a multimachine system 
in the form 

i = [A] y; + [B) u ; ^(o) = (1.18) 

is presented. The system model is transformed into the 
singular perturbation form as 


k 


^11 * 
i 

Ai2 


X 




j 

j. 

^21 1 



z 

+ 

j 



with x(o) = and z(o) = ^ , by properly identifying the 

perturbation parameter in the system equations. A logical 

T T T 

criterion for grouping the state variables as [x , ^ ], 
representing respectively the slow and fast variables, is 
presented. The technique first reduces the model (1.19) by 
neglecting the fast phenomena and then accounts for the effect 
of the fast phenomena by the addition of a ’boundary layer 
correction'. The decomposition of (1.13) into the fast and 
slow subsystems helps in understanding the structural proper- 
ties of the system and gives computational advantages as 
compared to the integration of (1.18). The computational 
benefits are derived mainly due to larger time step size 
permissible in the numerical integration of the slow and 
fast subsystems. A single machine-infinite bus and a three 
machine example are presented to illustrate the application 
of the theory to dynamic stability studies of power systems. 
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The singular perturbation theory can, in general, be 

extended to nonlinear system models. In Ohapter 5, an attempt 

is made to solve the transient stability problem using singular 

perturbation technique. The basic requirement in the singular 

T 

perturbation theory is to separate the variables ^ in the model 
T T 

as [x ,£ ] as the slow and fast state variables and z should be 
expressed as a function of x and u. This imposes restrictions 
on the selection of z variables in a nonlinear model of power 
system. Itirthermore , the computation of the improved response 
involves the solution of sets of differential equations having 
variab-le coefficients. A single machine-infinite bus 
power system example is solved illustrating the features of 
the technique. 

The concluding and the sixth chapter highlights the 
ma,-ior contributions of the thesis. Certain problems encountered 
during the course of this thesis are discussed and the suggestions 
for future work are then given. 



CHA.PTER 2 


IDETTTIPIOATIOTT OP GOHERMOY 

2.1 IRTRODUGTIOR 

A critical review of the existing methods for coherency 
analysis has been presented in Section 1.3.2 of the preceding 
chapter. Towards the end of the review a need for an improved 
criterion for identifying coherent generators has been brought 
out. In this chapter, an attempt is made to develop such a 
criterion and validate it on practical poAA/er systems. Before 
we describe the proposed technique, we define coherency between 
generators and coherency measures. 

a) Maximum Angular Excursion Measure [11] : For all practical 
purposes, two generators, say ’i' and *1’,are defined as 
coherent if 

1A6. (t) - A6.(t)| < e for o < t £T (2.1a) 

where = rotor angle of generator 'i' 

T = total transient period 
e = tolerancei 

The value of ^ is usually taken as 2,5° or 5°. 

h) Root Mean Square Measure [11]: For all practical purposes, 
two generators 'i' and 'j’ defined as coherent if 



where U is the tolerance. 


(2.1b) 
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The maximum angular excursion measure accounts for 
the maximum deviation in the rotor angles between the machines 
during the transient period. The r.m.s. coherency measure 
reflects the change in the transfer of energy between the 
generators during the transient period. The former measure 
is generally preferred in the formation of groups t11l because 
it ensures that there is no large difference between the rotor 
angle of the equivalent generator in the reduced system and 
the corresponding angles of the individual generators belonging 
to the coherent group in the full scale system. 

The desirable characteristics of the method for identi- 
fying coherent generators are that (i) the method must be 
simple to apply and efficient in terms of computer time, 

(ii) the results must be reliable, (iii) human judgement must 
be minimised and (iv) the method must be adaptable for multiple 
fault locations. 

In the proposed method, a new concept of coherency 
index is introduced and used as a criterion for identifying 
the groups of coherent generators. The linearised swing 
equations of the system are cast into the usual state space 
form, with the forcing function reflecting the fault under 
consideration [241. Then a closed form solution of the system 
equations is obtained in terms of the decoupled modes, k 
coherency index between two machines , which reflects a bound on 
the maximum angular excursion between the machines during the 
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transient period, is derived from the above closed form 
solution and used as a criterion for coherency analysis, 

A systematic procedure for grouping the generators is 
presented. The method does not require the storage and 
comparison of sv/ing curves. The technique is easily adapt- 
able for multiple fault locations i.e. for changes in the 
fault location in a study system or for changes in the geo- 
graphical configuration of the study system itself. The 
proposed technique is validated on typical power system 
examples. 

2.2 ASSUMPTIOITS 

The following assumptions, which are same as those in 
[241 , are made in dev-elopirtg the mathematical model for 
coherency analysis, 

(i) Groups of coherent generators can be determined by 
considering the linearised mathematical model of the power 
system. This is justified from the observation in the usual 
transient stability studies that the groups of coherent 
generators are independent of the size of disturbance. This 
assertion can be easily confirmed by considering a fault and 
observing that the coherent behaviour of the generators is 
not altered by changing the fault clearing time. 

(ii) Groups of coherent generators are independent of 
the detailed modelling of a generating unit. Therefore a 
classical synchronous machine model is considered. This 
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assumption is "based on the argument that though the control 
equipments have appreciable effect on the damping in the 
system, the natural frequencies, v^hioh predominantly deter- 
mine the coherent behaviour of generators, do not alter 
appreciably [73]. 

(iii) The response of a faulted system is approximated 
by the response of the unfaulted system by increasing the 
mechanical input power of each generator by an amount . 
(accelerating power at t = o'*’) for the time equal to fault 
clearing time only. 

For typical fault clearing times and inertia constants, 
the accelerating power of each generator in the faulted period 
remains nearly constant. An approximate variation of the 
accelerating powers ;in the post-fault period can be reproduced 
by applying a step input to each generator as shown in 
Pig. 2.1. Other implicit assumption in the above modelling 
of an electrical fault is that the effect of changes in the 
network: configuration of the study system and of the subse- 
quent changes in the load flow on the coherent behaviour of 
the machines fs ignored. Such a modelling of an electrical 
fault is valid for small fault clearing times and large 
systems. 

In hypothesis (iii), accelerating power of each 
generator at t=o’^ is computed from the voltage profile at 
that instant. The required voltage profile in turn is 




Time in seconds 


Fig.2*lQ Approximate variation of aoceterating 
power of a generator in faults 
system 





"s— Electrical Pow^ 



' te ' ■ ■ 

Time in seconds 

Fig. 2*1 b Equi volent reproduction of CKxrelerating 
power in the unfaulted system 
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computed either hy a load flow solution at t = o"^ or by a 
simpler noniterative procedure given below. 

(i) The constant voltage behind reactance together with 
the reactance of each generator is converted into a 'N'orton’ s 
equivalent. The magnitude of the current source is constant. 

(ii) Treating the transient reactances now as shunt 
branches, the bus admittance matrixlY] is computed. 

(iii) The network equations I =[Y]y are then solved for 
voltages by simulating the fault and the new generator powers 
are computed [12]. 

2.3 MATHEMA-TIOIL MODEL OP POWER SYSTEM 


In a power system having ’n’ number of generators, 

swing equation of the ith generator is given by 
,2 . 


Ml 


+ D, 


dd. 


dt‘ 


i dt 


n 

mi 13 

1 = 1 . 2 , 




. ,n 


( 2 . 2 ) 


Linearizing (2.2) around an operating point, we have 




d^( 


dt 


A6^) 
2 


-f D4 


d(A6^) 
dt' 


= Ap 


mi 


n 


y E.E.(-G. . 
3tl " 


Sin 


6? . 

ID 


1 Bi3 Cos (2.5) 

i = 1 , 2, . . . ,n 

Under the hypotheses (i) - (iii) and considering the nth 
machine as reference generator, equation (2.3) is cast into 
the state space form in the usual manner for the faulted 
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; t > tg 


with X(t ) = X (2.5) 

*** “ •"tJ 


and post-fault state of the system as follows. 

X = [Al X + iBl u ; o< t <tg with X(o) = o (2.4) 

X = [Al X 

where X = state of the system at time ’t^’ evaluated from(2.4) 

accelerating powers. 

The vector X and the matrices [A] and iBl are defined for the 
cases of uniform and non— uniform damping separately [yi ^ • 

Uniform Damping : Bor 0 = = B 2 /M 2 = ... = X, [A] 

and [Blare defined as 

iT 


• • • > 


Am . 1, 

n-1 ,n ' 


X= A 62 ^, •••» ^^n-1,nl^'^1n»^'"2n’ 

(2n-2) dimensional vector 


[A] = 


n-1 


n-1 


n-1 
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The elements of [ a’ ] are given by 




- I, rlf 

E. 


n 1 

\ 


(-G .Sin 6 - + B . Cos 6 . ) 
^ ni ni ni ni^ 


^id 


IS ]E) EE 

4^(-G,.^Sln6?, + B,^Cos4?,) - -^(-S_,Sin6°, +B„,Cos«°J 




1C 


1C 


1C 


1C 


M 


n 


nc 


nc nc 


nc 
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The primary motivation for selecting different state 
variables in the above two cases is to avoid the occurrence 
of a zero eigenvalue. The selection of state variables is 
also consistent with the concept of minimal realization as 
explained in Ref. [711. In the subsequent analysis, we 
consider the case of uniform damping only for the sake of 
simplicity. However, the analysis is valid even for non- 
uniform damping by taking appropriate [A] and [B] matrices. 

The elements of [A’ ] depend upon the system parameters 
and the pre-fault operating point. The expressions for the 
elements of [Aland [Bl are general. However, if as in 
Ref. [241 , the linearisation is done around the no.-load 
operating point i.e. the flat voltage profile, then the 
elements of [a’I can be directly derived from Stott^ s [B*] 
matrix [ 25 l and the inertia constants. 

In equations (2,4) and (2.5), the elements of [Al , 
[Bland the state vector X depend upon the choice of the 
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reference generator, whereas the vector u, which reflects the 
fault Under consideration, does not depend on it. For the 
sake of convenience and in order to establish a consistent 
relationship between X, [A] and [B] with u, we will always 
select a generator close to the fault location as the reference 
generator. Such a choice of reference generator gives the 
state variables as the relative rotor angles and velocities 
with respect to a generator, which is sensitive to the fault. 
The selection of reference generator near a fault location 
also forms the basis of the decomposition of a power system 
described in the next chapter. 


2.4 ANALYTIOAI SOLUTION 

If the eigenvalues of matrix [ A ] are distinct (which 
is generally the case with practical systems) the corresponding 
eigenvectors are linearly independent and the modes of oscilla- 
tions are decoupled from each other. Let the matrix tAl have 
’m’ pairs of complex and (2n-2m-2) real eigenvalues. The 
response of the system in the faulted state i.e. the solution 
of equation (2.4) is given as (Refer Appendix A for details) 

X^(t) = [P] z^(t) (2.6) 

where Xf = state of the system in the faulted period 

2m columns (2n-2m~2) columns 


r H t 

[P] =[£-]» > £ 2 » £.2 ’ 


t If 


£2m+1 »^2m+2’ •* ’£2n-2^ 


(2n-2)x(2n-2) matrix 
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The eigenvectors and reciprocal basis vectors are 
normalised and satisfy the relationships [70] between them 
as given below. 



Ii> 

= 1 , 

T? 

il> 

= 1 

<xh 

r" > 
ii 

= 0 

ft 

<li» 


= 0 


Similarly the response of the system in the post-fault period 
i.e. the solution of equation (2.5) is given as (Refer 
Appendix ^ for details) . 


Xpf{n = tsi is(t) 


(2.S) 


where 


= state of the system in the post-fault period 


2m columns 


(2n-2m-2) columns 


[S] = [s^, , S2, ^2' •••' ^ ' — 2ra+1'* — 2m+2^' 

t 

(2n-2)x(2n-2) matrix. 


»~2n-2^ 



2m 

(2n-2)-diia-eQaiDn.al ve-otor 

(2n-2m-2) 
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For the modes corresponding to complex eigenvalues, i = 


f 




ff . 



For the modes corresponding to the real eigenvalues, 
i = 2m+1 , 2n-2, ~ ^ — i^ 

The scalars 5,^ and are defined as 


Si'. 

1 



h> 




X > 


(2.9) 


The first 2m columns of [P] and [S] correspond to 
the complex modes and the remaining (2n-2m-2) columns 
correspond to the nonoscillatory modes respectively. The 
scalars k:[ and Ic^' in (2.6) and and in (2.9) depend 
upon the fault under consideration. 

2.5 IDMTI FIXATION OF GOHEREUT GENERATORS 

From the above theory -we now summarise the possible 
ways of identifying coherent generators including the con- 
ventional methods of coherency analysis. 

(i) Base Case Swing Curves ; In this case, the swing curves 
are obtained by numerically integrating the base case swing 
equation (2.2), The network changes due to the initiation 
and clearance of a fault are considered. -These swing curves 
are compared and the groups are determined by using the 
criterion (2.1a), 

(ii) Linear Simulation: In this method, the linearised 


swing equations (2.3) are numerically integrated over the 
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transient period "by using the trapezoidal rule and a fault is 
reflected as given in hypothesis (iii). The resulting swing 
curves are compared. 

(iii) Use of Closed Form Solution : The swing curves for a 
given fault are computed "by substituting for different values 
of ’ t’ in (2.6) and (2.8) and compared. 

In the next section we present a direct method for 
coherency analysis which does not resort to comparison of 
the swing curves. All the above methods of coherency 
analysis are critically compared in Section 2.13. 

2.6 DIRECT METHOD OP COHERENCY ANALYSIS 

Por the typical small fault clearing times obtainable 
in power systems, the differences in the rotor angles at 
t=o’'' and t=tg are very small. Hence it is sufficient to 
examine only the response (2.3) in the post-fault for 
;4s certain ing the coherency among the generators. To classify 
the generators ’i’ and ' as coherent, the maximum angular 
deviation between the state variables and should 

be examined. Prom eq.uation (2.3), we have 

"" ^-i’ ^ ^ ^ (2.10) 

where r . and r , are the row vectors of matrix ISI correspon- 
ding to the state variables A6j^ and a^^^ respectively. The 
expression f or A 6 • ^ is then written from (2.10) as 


-<(r^ -rp, £g(t)> 


( 2 . 11 ) 
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TTow consider the theoretical coherency between the 
generators 'i’ and ’ 3 ' i.e. A 6 j_j(t) = o. The condition for 
theoretical coherency is satisfied if — 3 II " "where 

II .ll represents the Euclidean norm. Since perfect coherency 
is rarely achieved, a practical coherency criterion, where 
certain tolerance is admissible, is needed. This criterion 
should have theoretical basis and physical justification. 

The following discussion provides the justification. 

Suppose that the eigenvalues of [A] are all complex 
with zero real parts, then ||£,-|| s^nd ||r.|j represent the 
bounds on the maximum amplitude of state variables A and 
A 6 ^ respectively. This point is explained as follows. 

(i) The length of vector ^g(t) is constant at any instant 
of time. Therefore, in a Tm-dimensional space, the vector 
£g(t) is rotating and always constitutes the radius of a 
hypersphere, 

(ii) A 6 . and A 6 . will have the possible maximum ampli- 

^ ^ in 3 n 

tudes when z„(t) is collinear with r. and r. respectively. 
Since the length of vector £g(t) is constant, llr^ll and 
Hare the bounds on the amplitudes of state variables 
A 6 ^^ and A 6 ^^ respectively. 

The above bounds are also valid if the matrix [A] 
has eigenvalues with negative real parts. With the classical 
representation for generators, usually the eigenvalues of lA) 
have small or zero real parts. This method of argument leads 
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2.6.2 Coherency Index 'h-^’: 

^ J 

Generators 'i' and ’j* are classified as practically 
coherent if 



II 1 1 ^ 

Max llrj^ll 

k=1 ,2, . . . ,n-1 


( 2 . 12 ) 


where h^^ is the coherency index between the generators ’i' 

and and ^e’ is the tolerance. In physical terms, the 

numerator represents how far the respective components r^^ 

and r^ differ from each other and thus gives a qualitative 

measure of deviation between and A6 . . The denominator 

in gn 

represents the normalising factor. 

a) Selection of e : In processing the swing curves for 
coherency, usually the tolerance of +2.5° or +5° is taken 
irrespective of the maximum angular excursion in the system. 
Since the coherency index is defined as a normalised measure, 
the value of ^ will depend upon the normalising factor 
Max||rjj.I|, k = 1,2,...,n-1. As a mle of thumb, the value of 
e is approximately selected by the following relationship 

e ~ Allowable tolerance in degrees 

Max II r,, II in degrees 

k=1,2,. ..,n-1 

2.7 FORMATION Of OOHERENI GROUPS 

The procedure for forming the groups of coherent 
generators uses a transitive process [12 ] i.e. if generator 
’a’ is coherent with generator 'b* and generator ’a’ is 
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is coherent "with generator ’ o* , then the generators ’h' and 
’ c’ are also coherent. In the grouping process, a comparison 
generator is defined for each coherent group and all other 
eligible generators are al-ways compared against this generator 
in order to determine whether they should be included in the 
same group or not. 

To start with all generators in the external system 
are marked as eligible and any generator (usually the 
farthest from the fault location), say ' i’ , is selected as 
the comparison generator of the group one. The coherency 
indices h. ^’s for all j's belonging to the set of eligible 
generators are computed. The first group of coherent 
generators is determined by applying the criterion for 
coherency (2.12). Then the generators included in the 
first group are deleted and the rest of the generators in 
the external system are marked eligible for forming the 
second group. The second group of coherent genei^tors is 
formed in a similar manner. In this way all possible groups 
of coherent generators in the external system are determined. 

2,8 COMPUTATION OP EIGENVALUES, EIGENVECTORS AND RECIPROCAL 
BASIS VECTORS 

In this section, we outline a method for computing 
the eigenvalues , eigenvectors and reciprocal basis vectors of 
[A] from those of the submatrix t A’ ], which is of dimension 
(n-1, n-1 ) (uniform damping). 



44 


Suppose- X is an eigenvalue of [ A 1 with the associated 

m m 

eigenvector , 2^2^' where and X 2 (n-1 ) -dimensional 

vectors. Then we have 

(2.13) 

This gives 

X2 = XX-i 

[a’] X^ = X(X+ 0) X-] (2.14) 

Hence, it follows that if X’ is an eigenvalue of [A’l , then 
there are two corresponding eigenvalues of [A] such that 


[ 0] 1 

[I ] 


2i . 



[ A’] 1 

L. i 

-ecu 


2^2 : 

= X 

: 

S2 


I 

X( X + G ) = X 

-0 +\/c^ -t- 4x’ 

or X = ( 2 . 15 ) 

2 

The following possihlities exist: 


X’ 

positive real 
negative real 

complex pair 


X 

positive real and negative real 

a complex conjugate pair or two 
negative real eigenvalues 

two complex pairs 


Our experience with Ca’ ] matrices for different systems 

shows that all the eigenvalues of [ a’ ] are usually negative 

2 

real. Moreover, since (0 +4X’ ) is generally less than zero 
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"because of small damping, "we will always have complex eigen- 
values of [Al. Then the eigenvalues, eigenvectors and reci- 
procal basis vectors of [A] can he easily derived from those of 
[A*] by using equations (2.15) and (2.14). The reciprocal 

basis vectors of [a'] are derived from the corresponding eigen- 
1 T 

vectors of [A ] . This facilitates the computation of 
coherency indices directly from the eigenvalues, eigen- 

X J 

vectors and reciprocal basis vectors of [A*] itself. 

The above technique is, however, not applicable if 
nonuniform damping in the system is considered. 

2.9 NUMERICAL EXAMPLE 

Before considering the problem of identifying the 
coherent gnerators for multiple fault locations, we now 
illustrate the above theory with the help of a typical power 
system in India. The power system belongs to the Uttar 
Pradesh State Electricity Board and consists of 71 buses, 

94 transmission lines and 13 generators. A single line 
diagram of the system is shown in Eig, 2.2 and the complete 
system data is given in Appendix B. The study system, which 
constitutes the western part of the system and the external 
system are defined as follows (Eig. 2.2), 

Study system = generators 1-5 
External system = generators 6-13. 

A three phase fault on bus 8 cleared in 0.1 sec. is consi- 
dered. The generator 4 which is close to the fault location 
is taken as the reference generator. 




machine sy 





47 


The matrices [A’j and [B] for the above system are 
given in Tables 2.1 and 2.2 respectively. The eigenvalues of 
[A'] , "which are all real and negative, are given below. 

-200.253, -3S.031 -116.525 -63.377 -140.665 -73.662 

-90.955 -117.436 -105.392 -110.137 -114.532 -109.394 

For the sake of brevity, the eigenvectors and the reciprocal 
basis vectors of [A'] are not given here. The vector u which 
reflects the fault under consideration is given by 

u^=[ 1 .2833, 1 .0313, 0.2153, 1 .4350, 0.4536, 0.5103, 0.6374, 
0.0794, 0.3455, 0.9142, 0.5072, 0.1544, 0.1575 ] 

The coherency indices are computed by using the procedure in 
Section 2.3. These coherency indices are given in Table 2.3 
and the table is written in accordance with the procedure for 
formation of coherent groups as explained in Section 2.7. 

The coherent groups with a tolerance e= 0.35 corres- 
ponding to 5® in the base case s"tudies are given below. 

Group 1 : generators 9, 11, 12 and 13 

Group 2 : generators 6 and 7 

These results are verified by plotting the swing curves 
from both (i) the base case simulation and (ii) the closed 
form solution given in equations (2.6) and (2.3). These swing 
curves for the generators in the above two groups are given in 
Pigs. 2.3 to 2.6. The base case swing curves indicate that 
the results obtained from the direct method are satisfactory. 
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Table 2.2 ; [BJ Matrix 
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Table 2.3 ; Coherency Indices 


Comparison 
Generator of 
Group 

Generators in 
External System 

12 

10 

8 

7 

6 

0.8261 

0.7119 

0.7831 

0.2616* 

7 

0.7879 

0.6554 

0.7552 

0.0* 

8 

0.5157 

0.6897 

0^0 

- 

9 

0.1643* 

- 

- 

“ 

10 

0.6507 

0.0 

- 

- 

1 1 

0.1 332* 

- 

- 

- 

1 2 

0.0* 

- 

- 

- 

13 

0.3379* 

- 

- 

- 


Maxllr^ll , k = 1,..,n-1 = O.2425 

e = 0.35 




%g. 2*3 Base case swing curves of generators 9,11.12 and 13 



6*06h 


52 



0-0 0-15 0-30 0-45 0-60 075 0-90 

Time in secs. ► 

Fig 2-4 Swing curves (closed form solution) of 
generators 9,11,12 and 13 



mu 


S4 



Fig. 2‘6 Swing c urves (closed form solution) of generators 
6 and 7 
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2.10 CHANGES IN THE EAULT LOCATION 

The prohlem of identifying coherent groups for different 
fault locations is an important problem in coherency analysis 
for the following reasons, 

(i) Consider a case where the concept of study system is 
not used and an equivalent is reconstructed for every individual 
fault location. In such case, whenever the fault location is 
changed, new groups of coherent generators are to be determined 
(discussion oa Ref, [24I), 

(ii) When the concept of study system and external system 
is used, we need a group of coherent generators valid for a 
set of fault locations in the study system for constructing 
the equivalent. In such a case, the groups of coherent genera- 
tors for individual fault locations in the study system are 
determined and then a set of coherent groups valid for all the 
fault locations is selected. 

(iii) Suppose the geographical configuration of the study 
system is altered and it is required to construct an equivalent 
for the new study system, then the step (ii) has to be repeated 
for new study system. 

Whenever a fault location is shifted, the following 
changes are incorporated in equations (2.6) and (2.8). 

(i) Accelerating power of each generator at t=o''' is 
recomputed resulting in a different forcing function u. 

(ii) Since a generator close to the fault is selected 
as the reference generator, a change of the fault location 
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may necessitate the change in the reference generator. The 
recomputation of nev? [A] matrix, its eigenvectors and reci- 
procal basis vectors is avoided hy the use of similarity 
transfomation described in the next section. The eigen- 
vectors and reciprocal basis vectors can be easily derived 
from the existing ones and a transformation matrix. 

(iii) These changes are correspondingly reflected in 
new values of and k? in (2.6). The new values of and 
are then computed from as given in (2.9). 

2.11 SIMILARITY TRAT^SIORMATIOM MATRIX 

Suppose we have already grouped the generators for 
a given fault with the reference generator, say ’ q.' . Mow 
it is required to group the generators for a different fault 
location and the corresponding reference generator is ’p’. 
Let the two state vectors be denoted by Xri — n 

two system matrices by ^ and [A^] respectively. For' a 
given pre-fault loading condition-^ the two matrices are 
related through a similarity transformation matrix. The 
transformation matrix is given by 

Xq = t T ^ 5p (2.16) 

where 



The components of ^ and X are defined in equation (2.4). 
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Then hy inspection [T] is obtained as 
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where [II is the identity matrix of appropriate dimension. 
Por p <q» matrix [T*] is obtained in a similar fashion. Then 

[ A 1 = [T]'"'' [A 1 [T] 

" , " ( 2.17 

[Bp] = [T]~' [Bq] 

The eigenvectors of matrix E-^p] can be derived from the 
eigenvectors of I the relation 
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2,12 ITOMERIGiLD EXAMPLES 

¥e now demonstrate the use of similarity transformation 
in the grouping of coherent generators for different fault 
locations. Consider the power system in Section 2.9. The 
coherent groups corresponding to two typical fault locations 


are eiven helow 
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a) Fault on Bus 15 (Same Study System): The study system is 
the same as in Section 2.9. The groups of coherent generators 
for a three phase fault on "bus 15 and cleared in 0.1 second 
are to "be identified. We have [A’] matrix, its eigenvectors 
and reciprocal "basis vectors with generator 4 as the reference 
generator. For the fault on tus 15, generator 5 is taken as 
the reference generator. The new eigenvectors and the 
reciprocal "basis vectors are derived by using the technique 
of similarity transformation given in Sections 2.10 and 2.11. 
The vector u of the accelerating powers of the generators at 
t = o’*” is given by 

u^ ={1.222, 0.902, 0,086, 0.105, 0.920, 0.514, 0.713, 

0.0S2, 0.554, 0,960, 0.516, a.157, 0.159 1 

The coherency indices are given in Table 2.4, The coherent 
groups for a tolerance 0.21 are as follows. 

Group 1 : generators 9, 11, 12 and 15 

Group 2 ; generators 6 and 7 

As indicated in Section 2.6.2 earlier, tolerance e is small 
in this case because of the higher value of MaxUr^jl , 
k = 1,2,..., n-1. 

By moving the fault location in the study system, we 
have ascertained that the above groups are satisfactory for 
all fault locations in the study system. 
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Table 2,4 ; Ooherency Indices 


Oomparison 
Generator of 
Group 

Generators in 
External System 

12 

10 

8 

7 

6 

6 

0.4025 

0.3870 

0.3749 

0.1732* 

7 

0.3770 

0.36t5 

0.3588 

0.0* 

8 

0.3881 

0.2536 

0.0 

- 

9 

0.1134* 

- 

- 

- 

10 

0.3027 

0.0 

- 

- 

11 

0.03311* 

0.0 

- 

- 

12 

0.0* 

- 

- 

- 

13 

0.1259* 

- 

- 



Max II 

II 

1,2,.. ,n-1 

= 0.399 


e = 0.21 
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on Bus 68 (Different Study System): The new study 
system is taken as 

Study system ; generators 8, 9, 11, 12, 13 
External system: generators 1-7 and 10. 

A three phase fault on bus 68 cleared in 0.1 sec. is consi- 
dered. The generator 13 is taken as the reference generator. 
ITew eigenvectors and the reciprocal basis vectors are computed 
using the similarity transformation given in Sections 2.10 and 
2.11. The vector u for the fault is given by 

u^ = t1.104, 0.779, 0,060, 0.023, 0.469, 0 . 521 , 0.724, 

0.082, 0.375, 0.997, 0.64S, 0.201, 0.4451 

The coherency indices are given in Table 2.5 and the coherent 
groups corresponding to e = 0.3 are 

Group 1 : generators 1 and 2 

Group 2 : generators 6 and 7 

c) Other Ifumerieal Examples: 

In order to demonstrate the applicability of the 
direct method of coherency analysis for large systems, we 
have tested the method on a power system given in Ref . [31 
consisting of 31 generators, 133 buses and 253 transmission 
lines. The coherency indices and the groups for different 
fault locations are reported in Ref. [741. The coherent 
groups agree with the base case swing curves. The coherent 
groups in this system corresponding to one fault location are 
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Table 2.5 : Ooherency Indices 


Oomparison 

Generator of 1 

Gronp 

Generators in 

External System 

3 

4 

5 

6 

1 

0.0* 

- 

- 

- 

- 

2 

0.1305* 

- 

- 


- 

3 

0.4227 

0.0 

- 

- 

- 

4 

0.7166 

0,4284 

0.0 

- 

- 

5 

0.5242 

0.6135 

0.3011 

0.0 

- 

6 

0.5526 

0.6260 

0.3039 

0.7033 

0.0* 

7 

0.5218 

0.6014 

0.7790 

0.6693 

0.2659* 

10 

0.3982 

0.6040 

0.7358 

0.6089 

. 0.6605 



Max j| r^ll 

, k: = 1, 

2j***yii[-** 

1 = 0.2747 


e =0.50 
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given in the next chapter, where a modified method for grouping 
of coherent generators is presented. 

2 . 1.3 DISCUSSION 

Some of the important aspects of the techniques for 
coherency analysis presented in this chapter are given helow. 

(i) The primary motivation for the analytical approach 
using the concept of modes was to eliminate the repe+i+ive 
computation of swing curves, when the groups valid for a 
set of faults in the study system are determined. This 
ohjective has been achieved hy the appropriate use of the 
similarity transformation (Sections 2.10 and 2.11), To illus- 
trate this point, the computation times for different methods 
to identify the groups in UPSEB system corresponding to three 
different fault locations are summarized in Table 2.6. It is 
evident from the table that the computation time for the base 
case swing curves and linear simulation is proportional to 
the number of fault locations. But in the direct method 

and in the closed form solution method the additional computer 
time requirement for a subsequent fault location is about 29.59^ 
and 59.69^respectively . 

(ii) The linearised swing curves exibit reduced amplitudes 
and increased frequency of oscillations as compared to the base 
case swing curves. This observation agpees with that in '[ 24 ], 

(iii) In the mathematical modelling of power system 
(Section 2.3), we have considered only electrical faults as 
the disturbances. It is possible to obtain a general 





Table 2.6 : Oomputation Times for Ooherenoy Analysis in 
UPSEB System. 


!TOt of Base Ca 

fault S-wing 

locations Ourves 
considered 


14.0 


Particulars 

Proposed 
Swing Curves 
from 

Closed Pom 
solution 

Teohniques 
The Direct 
Method 

linear 

Simula- 

tion 

(a) 

4.05 

4.05 


(1) 

0.68 . 

0.68 


(c) 

1 .90 

1 .90 


(d) 

1.16 

- 


(e) 

- 

0.04 


Total time 

7.79 

6.67 

6.82 

(includes 
1 ,9 sec. 
_for_(c)_ 

(f) 

0.03 

0.03 


(c) 

1 .90 

1 .90 


(d) 

1.16 

- 


(e) 

- 

0.04 


Total time 
(inclusive 
of the 
above) 

10.38 

8.64 

13.64 

(f) 

0.03 

0.03 


(o) 

1.90 

1.90 


(d) 

1.16 

- 


(e) 

- 

0.04 


Total time 

13.97 

10.61 

20.46 


28.0 


42.0 


(inclusive 
of the above 
two) 


Oon tinned . . 
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ITotes: 

(l) All entries in the table are GPU times in seconds on 
Dec System 10. 

(ii) Swing carves in the base case and the linear simulation 
method are plotted for 1 sec. only with an interval of 
0.01 and 0.05 sec. respectively. 

(iii) Gomputer time for comparison of swing curves/coherency 
indices which is usually very small, is not included. 

(iv) Gomputer programs for the base case transient stability 
and network reduction in step (a) do not use sparsity 
technique whereas the programs for computation of the 
accelerating powers in step (c) and the linear simulation 
method make use of sparsity technique. 

Particulars 

(a) Network reduction and computation of [A’] matrix. 

(b) Computation of eigenvectors and reciprocal basis 
vectors of t A' ] . 

(c) Computation of accelerating powers of the generators 
at t = o+. 

(d) Gomputation of swing curves from the closed form 
solution . 

(e) Computation of coherency indices. 

(f) Computation of the new eigenvectors from the existing 
ones and a similarity transformation matrix. 



66 


mathematical model of the system, -which also considers other 
types of disturbances such as loss of generation, tripping of 
a transmission line etc. In fact such a state space model is 
given in Ref. [27]. The theory presented in this chapter 
applies to such a generalised model as well. 

(iv) In the grouping of coherent generators (Section 2.7), 
the selection of comparison generator of a coherent group is 
arbitrary and the process is transitive. Therefore, for a 
given fault, the groups of coherent generators need not be 
unique. However, the size of the equivaliaed external system 
will not be affected. 

(v) The direct method gives a more physical insight into 
the dynamics of the system than the procedure involving 
processing of a large number of swing curves for coherency. 

(vi) A possible criticism of the method may be computation 
of eigenvectors and reciprocal basis vectors for large scale 
systems. Although not explored in this thesis, the concept 

of dominant modes may alleviate this difficul-ty, 

2.14 OOHOLUSION 

In this chapter we have presented two techniques for 
identifying coherency (i) by comparing the swing curves 
obtained from closed form solution and (ii) the direct method 
using coherency indices. The concept of similarity transfor- 
mation has been used to eliminate the repetitive computations 
involved in determining coherency for different fault locations 
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In the study system. This fact makes the method attractive 
as compared to existing methods. Validity of the method has 
heen tested on realistic power systems. 



CHAPTER 5 


DEGOMPOSITIOH OP POWER SYSTEMS AFP GOHSTRUCTIOH OP 
GOHEREHGY BASED EQUIVALENTS 

3.1 IHTROPUGTIOH 

Another prohlem, in addition to dynamic equivalence, 
facing power engineers is how accurate the model for a 
generating unit should he as one moves away from the fault 
location. This can he considered either as a step prior 
to dynamic equivalencing or a problem in itself. Por a 
long time, it has been the common practice to represent 
the generators- ’away* from the fault location with leaser 
detail than those in the ’vicinity* of it. Such a classi- 
fication of generators is not precise and is based upon the 
experience with the system studies. A logical basis for 
determining the levels of complexity in modelling the 
generators is the electromechanical interaction between the 
machines during the transient period. The electromechanical 
interaction between -two machines can be judged by (i) the 
maximum angular excursion or (ii) the synchronising power 
flow between the machines. In the litera-ture, the admittance 
and reflection distance measures [ 22 ] have been proposed for 
the purpose of modelling the generators. While the admittance 
distance is purely a static measure , the reflection distance 
does not reflect the electromechanical interaction between 
the machines during the entire transient period. 




69 


In this chapter, we (i) present a new method for the 
decomposition of a power system for the purpose of modelling 
the generators, which incorporates the electromechanical 
interaction between the machines during the transient 
period, (ii) establish a link; between the decomposition 
of power system and the direct method of coherency analysis 
presented in Section 2.6 for grouping of coherent generators 
and (iii) present a case study for a realistic power system 
to illustrate the principle of equivalen cing. The relevant 
equivalising procedure is also described briefly. 


3-2 ELETTROMEGHAWIOAL DISMOE (MD) 

We define MD between a generator in the external 
system and the reference generator in the study system as a 
on the maximum angular excursion between the two during the 
transient period. In terms of the mathematical development 
of the previous chapter, the MD between the generators ’i’ 
and ’n’ is defined from equation (2.10) as 


^in 


lUill 

Max. 


k;=1,2,.. .,n-1 


(5.1 


where generator ’i' belongs to the set of generators in the 
external system and ’n’ is the reference generator in the 
study system. as defined in equation (3.1) gives a 

per unit bound on the maximum angular deviation between the 
generators ’i' and 'n*. 


bound 
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3.3 DECOMPOSITION OP POWER SYSTM 

Based on the concept of EMD, a power system is 
divided into different regions for a given fault location. 

The division is done with respect to a reference generator, 
which, as stated in Section 2.3.1, is selected close to the 
fault location. The study system, by definition, constitutes 
region I, which corresponds to inner circle defined in 
Ref. [22], The proposed criteria for the demarcation of a 
power system into different regions is the following: 

Region I : Study system 


Region 

II 

‘ 

d. < 
in 

b2 

Region 

III 

: , b2< 

< 

in — 

^5 

Region 

IV 

: b^ ^ 

^in^ 

1 


where , ^2 and are the known constants, which define 
the boundaries between the regions. 

Since the transient behaviour of the generators in 
a study system subjected to the local disturbances is to 
be investigated in detail, the level of complexity in 
representing the dynamics of generators in the study system 
should be the highest. The generators in region IV being far 
away (in the sense of EMD) from the reference generator can 
be represented by either classical models or the dynamics 
of the generators in this region may be altogether neglected. 
The degree of details in modelling the generators in regions II 
and III is suitably graded in between these two extremes. 
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3.4 SEIEOTIO¥ OE THE BOUTEDARIES 

After decomposition of power system into regions, 
coherency analysis follows. Since the generators belonging 
to different regions are modelled with different degrees of 
details and difficult to aggregate [12], a coherent group 
spread over two or more regions is undesirable. Therefore, 
while selecting the boundaries !£ and b^, we have to ensure 
that a generator in one region will not be coherent with the 
generator in another region. This is done as follows. 


Coherency between the generators ’i’ and ’j’ implies 

that 



(3.2) 


i.e. the generators are equidistant (in the sense of EME) 
from the reference generator. However, the reverse may not 
be true. Hence (3.2) forms a necessary condition for 
coherency because . in the definition of EME, the sign of 
the deviation is not considered. Thus if and 

do not s tisfy (3.2), it is uaranteed that the generators 
' i’ and ’j' will not be coherent. The adjustment of 
boundaries is done so that the generators satisfying (3.2) 
belong to one region only, first the generators in the 
external system are arranged in the order of increasing 
EME’s, When a boundaiy b^ is selected, we will verify that 
the two generators on opposite sides of the boundaries and 
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close (in the sense of EMD) to it are not coherent hy comput- 
ing the corresponding coherency index given in equation (2.12). 
The selection of values of h^’s is, however, left to the system 
planner. Such a division of powar system facilitates the 
grouping of coherent generators in each region separately. We 
now illustrate the above technique of decomposition of power 
system and the subsequent grouping of coherent generators with 
two numerical examples, 

3.5 NUMERICAL EXAMPLES 
3.5,1 Exampl e 1 ; 

Consider the power system and the fault location 
given in Section 2.9. Generator 4 is taken as the reference 
generator and the EMD’ s of the generators in the external 
system are given in Table 3.1. 

Table 3.1 : Electromechanical Distances 


Generators 
in the 

External System 

lln 

Remarks 

7 

6 

0.6873 

0.7052 

Region II 

10 

0.7486 

^10,13 ^ 0-5700 

13 

0.8338 

e = 0.35 

8 

11 

0.9411 

0.9605 

Region III 

12 

9 

0.9863 - 

1 .0 



MaxllrJI = 0.2425 




The power system is divided into three different regions as 
shown in Fig.5-.1 using the techniq,ue given in the previous 
section. 
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Region I : Study system, generators 1-5 
Region II ; generators 6,7 and 10 

(0.65 < d^^< 0.75) 

Region III: generators 3,9,11,12,13 

(0.75 < 1 ) 

The coherent generators in each region are grouped separately 
using (2.13) as given in Tables 3.2 and 3-3. The coherent 
groups are as follows: 

Region III : Coherent group 1 : generators 9,11,12,13 
Region II : Coherent group 2 ; generators 6 and 7 

Table 3.2; Grouping of Coherent Generators in Region III. 

Comparison generator 

of group 12 8 

Generators 
in Region III 


8 

0.5157 

0.0 

9 

0.1643* 



11 

0.1332* 

— 

e = 0.35 

12 

0.0* 

- 

Maxllr^il = 0.2425 

13 

0.3379* 

— 

k=1 ,2, . . . ,n-1 
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Talkie 3,3 : G-rouping of Coherent Generators in Region II 


Comparison generator 
of group 10 

6 

Generators in 

Region II 



6 

0.7119 

0.2616* 

e = 0.35 

7 

0.6554 

0.0* 

Maxjjr^j] = 0.2425 

10 

0.0 

- 

*1 • • • ytl** 


3.5.2 Example 2: 

In this example, the decomposition and the subsequent 
coherency analysis of a pover system consisting 31 generators, 
132 buses and 258 transmission lines is presented. The single 
line diagram of the system is given in Fig. 3.2 and the system 
data is taken from Ref.[S], which for the sake of brevity is 
not reproduced here. The following case is considered. 

Study system : generators 4, 5, 7, IS, 19 and 20 

External system : rest of the generators 

Pault ■ : a three phase fault on the terminals 

of generator 5 cleared in 0.1 sec. 

Reference : generator 5. 

Decomposition t The electromechanical distances for the 
above case are listed in Table 3.4. The proposed decomposi- 
tion of the system is given below and is shown in Fig, 3. 3. 

Region I : Study system, generators 4,5,7,13,19 and 20 





IT 


Table 3.4 : Electromechanical Distances 


Generators in 
External System 

d. 

in 

Remarks 

28 

0.1290 


21 

0.1290 


S 

0.1301 


10 

0.1327 


22 

0.1331 


31 

0.1352 

Region II 

30 

0.1469 


11 

0.1499 


12 

0.1542 


23 

0.1840 


27 

0.1973 


26 

0.2060 

^26,29 0*2846 

29 

0.2165 

e = 0.12 

6 

0.2234 


9 

0.2244 


25 

0.2306 


1 

0.2351 

Region III 

13 

0.2405 


14 

0.2454 


24 

0.2703 


16 

0 . 2880 


15 

0.2977 

^ = 0,6551 

3 

0.5665 

e = 0.12 

17 

0.6565 

Region IV 

2 

1.0 



MaxJIr^j.|| = 0.8583 
k=1 ,2, . . . ,n-1 






group 

F = Fault location 

Fia 3*3 Decomposition of 31 machine system 
^ and formation of coherent groups 


79 


Region II : generators 28, 21, 8, 10, 22, 31, 30, 11, 
0.1 < 0,21 12, 23, 27 and 26 


Region III ; 
0.21 < d. <0.3 

II 


generators 29, 6, 
and 15 


9, 25^, 


13, 14, 24, 16, 


Region 17 : generators 3, 17 and 2 

0.5<di„<1.0 

We have ensured that there is no coherent group spread over 
two regions by computing the coherency index for two gener- 
ators close to each boundary as shown in Table 3.4. 

Groups of Coherent Generators ; The coherency indices for. 
the generators in the regions II and III are given in 
Tables 3.5 and 3.6 respectively. The groups of coherent 
generators for the tolerance e= 0.12 corresponding to +5° 
in base case studies are given below and shown in Rig. 3.3, 

Region II group 1 : generators 30, 28, 22, 21, 10 and 8 
group 2 : generators 31, 12, and 11 

group 3 : generators 26 and 27 

Region III group 4 : generators 29 and 25 

group 5 : generators 24 and 13 

In region 17, there is no coherent group as is evident 
from the MB's in that region. The above coherent groups are 
verified from both (i) the base case swing curves and (ii) the 
swing curves from the closed form solution. As an example. 
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Tatle 3.5 i formation of Coherent Groups in Region II 


Comparison generator 
of group 

Generators in 
the region II 

30 

31 

26 

23 

8 


0,0963* 

— 


— 

10 


0.0678* 

- 

- 

- 

11 


0.1746 

0.1014* 

- 

- 

12 


0.1864 

0.1135* 

- 

- 

21 


0.0964* 

- 

- 


22 


0.0940* 

- 

- 

- 

23 


0.2162 

0.1894 

0.2349 

0,0 

26 


0.2249 

0.1855 

0.0* 

- 

27 


0.1977 

0.1739 

0.0682* 

- 

28 


0 . 1125 * 

- 

- 

- 

30 


0.0* 

- 

- 

- 

31 


0.1516 

0.0* 

- 

- 



Max II 

EkII = 

0.8583 



k; = 1,2,... ,n-1 
e = 0.12 
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Tatile 3.6 : Pormation of Coherent Groups in Region III 


Comparison 
generator 
of group 

Generators 
in Region 
III 

25 

24 

16 

15 

14 

9 6 

1 

1 

0,3602 

0.3048 

0.4477 

0.4691 

0.2256 

0.2747 0.3591 

0.0 

6 

0.2413 

0,4532 

0.2380' 

0.2125 

0.3192 

0.8907 0.0 

- 

9 

0.3453 

0.3553 

0.4370 

0.4573 

0.2436 

0.0 

- 

13 

0.4201 

0.1058* - 


- 


- 

14 

0.3169 

0.3100 

0.4190 

0.4395 

0.0 

- 

- 

15 

0.4297 

0.5431 

0.3406 

0.0 

- 


- 

16 

0 .6660 

0.6790 

0.0 

- 

- 

^ - 

- 

24 

0.4580 

0.0* 

- 

- 

- 

- 

- 

25 

0.0* 

- 

- 

- 

- 

- 

- 

29 

0.1200* - 

- 

- 

- 


- 


Max. II rj^ll = 0.8583 

k = 1,2,... ,n-'1 

e = 0.12 
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the "base case and the closed fom solution swing curves for 
the generators in coherent group 1 are given in Pigs, 3.4 
and 3.5 respectively. 


3.6 OOTTSTRUOTIOF OP COHERENCY BASED EQUIVALENT 

After classifying the generators in the external 
system into different coherent groups, the next step is to 
replace each coherent group hy an equivalent unit. In this 
section we describe briefly the equivalising procedure of 
Systems Control Inc. tl2,13 1 which forms the basis for the 
case study presented in Section 3.7. 

3.6,1 Reduction of Coherent Generator Buses: 

In this method, the terminal buses of the generators 
in a coherent group are connected together instead of combi' 
ning the internal buses of the generators as in[7,8]. Such 
an approach decouples the problem of network reduction from 
that of dynamic aggregation of the generating units. It 
also facilitates the equivalent to be compatible with the 
existing transient stability programs. 


The network equations of the system are arranged in the 


following manner. 



“Ia‘ 





ib 

= 

'VI 



Ic 


[01 


where ' 

c’ represents 

’q’ numb 


[01 

■ 

"v " 


—a 



ib 



V 




er of terminal buses of the 
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generators in a coherent group, ’b’ includes the boundary 
buses of the group, which are ’k’ in number, and the buses 
in rest of the network are included in *a*. Since the 
replacement of the buses in ’ c’ is not going to affect the 
structure or the parameters in ’a*, we wiH. consider only the 
buses ’b’ and ’ c' in the reduction process. The network 
equations in (3.3) for groups ’b’ and ’c' for the purpose 
of reduction are 


• 


• 

Ik 


T 

■ tk+1 

• 


• 

^k+q 



^11 

* 

•• ^Ik 

« 

1 ,k+1 

I . 

•* ^1,k+q 

• 

• 

^k1 

• 

•• ^kk 

I • 

j ^k,k+1 

• 

* * ^k,k+q 


\+1,1 *• ^k+1,k* ^k+1 ,k+r* ^k+1 ,k+q 


fir. 


hr. 


hr 


k+1 


hr 


k+i' 
(3.4) 


*^k+q,1 ** ^k+q,k I ^k+qjk+l*" ^k+q,k+q 

Since the last ’q’ buses in (3.4) are replaced by an equivalent 
bus ’ t’ , equation (3.4) becomes 

I 1 r 

Y. 


^11 

• 

••• ^ik j 

• 1 

^1t 

• 

• 

^kl 

• 1 

* * • ^kk 1 

• 

^kt 

^t1 

* * • ^tk ! 

^tt 


'1 


7, 


7 , 


(3.5) 


The correspondence between (3.4) and (3.5) is established by 
the principle that(i) the total power of a coherent group 
should be preserved in the equivalent and (ii) the power 
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flovs over the interconnecting lines should he preserved 
in the equivalent. 

The expressions for the unknown quantities in (3.5) 
are then determined as ( 12 ] 




^tb 

= 2. 
i=k+1 

k+q 

= 1 
i=k+1 

^ib 

^bt 

^ib -- 


k+q 

k+q 

^tt 

= I 

I 

i=k+1 

d=k+1 




(3.6) 


IlZI 

■'t 


where ’h’ is a boundary bus between the group and the rest 
of the network. The elements in (5.6) are readily determined 
if i.e. the voltage of the equivalent terminal bus is 
known. The choice of is arbitrary and is usually taken as 
the average of the all terminal voltages in the group. 

ICH-q 0^ 


k+q .... ^ 

IVtl = I ' e = r (5.7) 

^ i=k+1 ^ ^ i=k+1 ^ 


iv.i 


In general, ^ indicating the presence of a 

phase shifting transformer between a boundary bus »b’ and 
the terminal bus ’ t’ . 

3.6.2 Dynamic Aggregation of Coherent Units: 

The process of determining the parameters and the 
equivalent generation from the corresponding parameters 
and the generations of the individual generating units 
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is termed as the dynamic aggregation of the generating units. 
Assuming the classical models for the generators, an equiva- 
lent generator representing a coherent group is defined by 
the following quantities and parameters. 


Generation : 

^g(eq) 



^g,k:+2 "■ ••• 

^ ^g,k;+q 


Qg(eq) 

" ®g,k+1 

+ 

‘^g,k+2 ^ 

‘^g,k:+q 

Load at the 

^S-(eq) 

= P 

£,k+1 

-f 

^jl,k:+2 + • • • 

^fi,,k:+q 

terminal bus: 






'^fi'(eq) 



'^£,k;+2 

*^Jl,k+q 

Inertia : 

= 

eq 


\+q 


Lamping 

Coefficient; 

L 

eq 

\4-2 

\+q 


Generator 

Impedance : 

1/z’ 

' eq 


''/^k+2 + ••• + 



The procedure for replacing a coherent group, by an 
equivalent generating unit and its modified interconnections 
as described in Sections 3.6.2 and 3.6.1 respectively is to 
be repeated for each coherent group; Since we have represented 
the loads by the constant impedances, the load buses are 
reduced by using the Gaussian elimination formula, 

3.7 A CASE STUDY 

In this section, we illustrate the principle of 
coherency based equivalent with the help of a case study 
for a typical power system. The power system considered 



is described in detail in Section 2.9 and the study system 
and the external system are shown in Fig. 2.2. The coherent 
groups valid for a set of fault locations in the study system 
as determined in Sections 2.9 and 2.12 are the following. 

G-roup 1 : generators 9, 11, 12 and 13 

Group 2 : generators 6 and 7 

The procedure of equivalising for the coherent group 1 is 
shown in Fig. 3.6(a) and (b) and that for group 2 in 
Fig. 3.7 respectively. Thus ,the equivalent of the system 
consists of 9 generators against 13 in the full scale system. 

A three phase fault on bus 10 cleared in 0.1 sec 
by opening the line 10-7 is considered. The responses of 
generators 1 and 3 for this fault are given in Fig. 3.8 
and Fig. 3.9 respectively for both (i) the full scale system 
and (ii) the equivalent system. 

The responses of the generators in the equivalent of 
the system match closely with the responses in the full 
scale system. The computer time required for simulating 
the equivalent is about 65^ of the time for the full scale 
simulation. This excludes the time for constructing the 
equivalent. A significant reduction in size and computer 
time is possible in large scale systems where the number of 
coherent groups is large. For example, as reported in 
Ref. [12], a power system having 337 generators is reduced 
to 166 generators and simulated in about 50f» of the time 
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Fig 3- 6a Coherent group no.1 in the full scale system 
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Fig. 3-7 (a) Coherent group no.2 in the full scale system 
(b) Equivalent representation 
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required for full scale simulation. However, in the present 
case the advantage is realised hy using the same equivalent 
for different fault locations in the study system. 

3.S GOHOLUSIOH 

In this chapter, we have presented a method for 
decomposition of a power system into different regions 
for modelling the generators. The criterion for decompo-r 
sition incorporates the electromechanical effects between 
the machines during the transient period* The demarcation 
of a power system into different regions has also been used 
to simplify the procedure for grouping of coherent genera- 
tors. The grouping technique and the construction of a 
coherency based equivalent have been demonstrated on a 
power system example. 



CH4PTER 4 


DYRMIO SIMPLIFIOATIOR U5IITG SINGULA.R PEHPURBATIQIT THEOBY 

(RYRAMIO STABILITY STUDIES) 


4.1 INTRODUCTION 

A model of much interest to power engineers is of 
the form 

£ = [A] + [B] u with £(0) = (4.1) 

This model is extensively used in the dynamic stability- 
studies and in the design of control circuits. The recent 
trend to represent the dynamics of a generating unit in 
detail incorporates the effect of many small time constants. 
Ignoring completely the effect of these time constants in 
the simplified model gives rise to a less accurate represen- 
tation of the given system. Therefore the technique for 
simplifying (4.1) should consider the effect of such time 
constants to improve the accuray. A need for an effective 
technique for this purpose has been brought out in Section 
1.3.3, in which the existing techniques for simplifying 
(4.1) have been discussed. In the proposed technique we 
make such an attempt by using singular perturbation theory. 

This chapter is organised along the following lines. 

Firstly a systematic algorithm to formulate the 
model (4.1) for a multimachine system using the detailed 
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representation of the generating units is presented. Then 
the model is written in the singular perturbation form (1.19) 
by treating the response of (4.1) as the perturbation of the 
response of the simplified system (reduced system). This 
needs the appropriate separation of the fast and slow state 
variables. The validity of the simplified system is ascertain- 
ed by the condition given by Vidyasagar [60], The subsequent 
improvement in the response is obtained by deriving the first 
order solution to (4.1) by using the method of asymptotic 
expansions. The technique results in solving a set of fast 
and slow subsystems separately and where the step size of 
integration could be larger than for the original system (4.1), 
The resulting computational saving makes this method attractive 
for model simplification. Single and three machine examples 
illustrate the technique, 

4.2 POWER SYSTEM MODEL 

A power system model for dynamic stability studies 
consists of a linearised set of differential equations 
(representing the dynamics of generating units) and a set 
of algebraic equations (describing their interconnections). 
These equations are cast in the state space form (4.1) by 
©liminating the non— state variables. We consider a synchro- 
nous machine, an exciter-voltage regulator system and a 
turbine-governor system as the components of a generating 
unit. A synchronous machine is represented by Riaz's hybrid 
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parameter model with equivalent rotor voltages as state 

variables [62,64]. The voltages A eJ and AE” are chosen 

as state variables rather than Ae" and AE^ as defined 

in [62]. The rotor voltages and AE^ proportional 

to the flux linkages of the damper windings are associated with 

the time constants T^^^ and T^^. The rotor voltagesAE" andAE" 

a q 

are then expressed as 


ae" 

q 



= % * ae" 

= ae' . AE" 


(4.2) 


The standard models for an exciter-voltage regulator system 
and a turbine-governor system as given in Ref. [ 65-68 ] are 
used. The details of modelling a generating unit are given 
in Appendix 0. In the main text of the thesis, we will use 
only one symbolic state space equation for each component 
of a generating unit. 


4.2,1 Generating Unit Model: 

Consider a power system having ’n’ number of 
generating units. The state space representation of ith 
generating unit with its own d-q axes as the reference frame 
is given by 


Li 


A , 
mi 

*^mel 

mti 


L 


1 

a 

B 

J 

AI . 
-mi 

2 Sml 

Li 

=: 

0 . 

emi 

A . 

ei 

0 


^ei 


“el 

-h 

iei 2 

-P 

i 





•^ti 


“ti 


° Iti 


i=1j2,,,,,n ( 4 . 5 ) 



9S 


where and 

the state vectors representing the 
dynami cs of the ith synchronous machine, its exciter-voltage 
regulator system and the turh in e-governor system respectively. 
These are defined as follows: 

=tAE;^, A,»., AS., aE;., aE^^, AE".1 

£ei = ^621’ ye31’ ^e4i' ^eSl ! 

4i = t ytii' yt2i' ytji 1 


^imi — gi defined as 


-flii ‘ di qi ■* ' — gi refi^ moi 


The details of equation (4.3) are given in Appendix 0 
and the state variables and are identified on 
Pigs. 0.2 and C.3 of Appendix 0 respectively. 


4.2.2 Transformation to the Network Reference frame: 

The following transformation is applied to refer 
the machine currents and voltages to a common network referen- 
ce frame denoted by D-Q axes [641 (Pig.4.1) 


r 1 


— 





Gos 6^ Sin. <5^ 





1 

-Sin 5^ 0os6^i 


^Qi, 

^ — 


— 

i 


1,2,. 


(4.4) 
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On linearisation, equation (4.4) becomes 


■where 




^Di’ ^^Qi^ 


,n 


[S,l = 


(bp 6° 


in 6? 
1 


Sin 


-S-irr 6° Oo« 6° 


J. 

-1 


1 °. 

qi 

-1°. 

di 


(4.5) 


4.2.5 TTetwork Equations: 

The transmission network equations describing the 
interconnections of the generating units are given by 


1— * 


^ ’■DI 




• 1 

* ’-Dn ■ 


"’’Qn: 

1! 

where 



-^11 

®11 ^11 


^n1 


B 


-B 


'n1 


n1 n1 


"ln 

®1n "ln 


G -B 
nn nn 


B & 

nn 




AE, 


QI 


AE 


'3>n 


AE 


nn 

_j U. 


Qn 


(4.6) 



= AE^^ + AE^^ 

(4.7) 

A E^. 

Qi 

= . E' , . AB", 


The inverse 

transformation fo r 

voltages and currents is 

given by 




= ^ ^^i + -i^ *^i 

■.It 

i=1,2,...,n (4.8) 
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where 


A^'. = (AS 


ft 

Di' 




. = (AS^^, AE;.]^ 




,0 

'Di 


T 


4.2.4 System Matrix: 

The system matrix is formed ty arranging the 

equations (4-3) for 'n’ numter of generating units and 

eliminating the non-state variables Al^ = AI^ 1 

“in -ini ^ -iun 

by using equations (4.5), (4.6) and (4.8). The resulting 
system equation is given as 


^ = [Al y; + [bI u , with. y^(O) = ^ 


( 4 . 9 ) 


where 

..T 


T T !P T T G? T 

^g2’ ’ • • » •^gn^ ^ '^gi “ ^^i’ ^ei’ -^ti^ 


u 


T 


T 


T 




The rotor angles in (4.9) can be expressed with respect to 
a rotating d-q axes frame of a reference machine, thereby 
reducing the order of the system by one in accordance with 
the procedure given in Ref . [61 ,643 . 

Without loss of generality, let the dimensions of ^ 
and u be denoted by *m’ and ’r’ respectively. Wow we are 
concerned with the dynamic simplification and simulation of 
( 4 . 9 ) using the singular perturbation technique. Before this 
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problem is taken up, "we briefly review the singular pertur- 
bation theory as applied to the initial value problem in a 
linear time invariant system. 


4.3 SIT<rGUMR PERTURBATIOF THEORY 

Suppose the system (4.9) is written in the singular 
perturbation form as 





1 


- 


— 

X 

. — ^ 

■ = 

Ai ^ 

1 ^12 


X 

, + 




1 

> 

(V) 

1 

1 ^22 


z 




with x(0) = x° and ^(O) = , 


( 4 . 10 ) 


where x and z, have dimensions *m^^ and ’m2’ respectively such 
that m^+m2 = Assume that 

(i) the perturbation parameter ’e' is small and > 0 and 

(ii) the matrix ^-^22^ nonsingular and stable. 

The asymptotic solution of (4.10) is an additive function 
of time ’ t’ and a stretched variable x = t/e. We seek the 
solution of (4.10) in the form [50 ] 

x(t, e) = X(t, e) + e£(T, e) 

zit, e) = 2(t, e) + X, e) 

where (X, Z) is the outer expansion (i.e. the solution away 
from t = 0) and ( e2» boundary layer correction, 

which is significant only near t = 0. X, 2, £ and all 
have asymptotic expansions as e -»-o as 
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(S, Z. £>2.) -- ( ! Xj f Z I <P, I 4 c^) 

( 4 . 12 ) 

Under the specified hypotheses (i) and (ii), the boundary 
layer terms £ and £ in (4.11) will decay fast as e -^-o 
(or T-»- c°) S’Hd (x, z_) will converge to (X, Z) » The equations 
for different terms in (4.12) can he derived by substituting 
the expansions for (X, Z) and (x, z) in (4.10) separately 
and equating the terms corresponding to equal powers of e 
on both sides. 

The first term (X^, Z^) in the outer expansion is the 
solution of the reduced problem, which is given by setting 
e = 0 in (4 »10) . 

^ ^R^ - ^ (4.13a) 

= .- IA“2 t^22 ®2^ ^ (4.13b) 

where 

[Aj^]= [All ~ “^1 2 ^22 "^21 ^ 

[ = [ B-j “>A.-| 2 ^22 ®2 ^ 

is obtained by solving a set of differential equations 

-^2. = ^-^22^ ^0 ’ = 2° - 2^(0) (4.14) 

dT , 

and T = t/ e 

Then the zeroth order solution of gi"^^ by 

x(t) = XQ(t) + 0(e) 
z(t) = Z^(t) + ^^(t) + 0(e) 


(4.15) 
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^Further inipjc 0761116111 : in th6 approximation to (x, z) can te 
obtain6d by doriving tho first order solution of (4.10) as 
given bslow. 


x('t) = + cl X-iCt) + £o(t)] + O(e^) 

z(t) = Z^Ct) + i^(T) +e[ Z^(t) + £^(t) ] + 0 (e 2 ) 


( 4 . 16 ) 


This, ho-wever, involves the solution of following additional 
sets of equations for (X^ , Z^) and (£^, ) with appropriate 

initial conditions. 


dx 

= EA. 

12 ' 



(4, 

.17) 

^(~) 

= 0 

and 


) is 

given by (4.14). 



« 




-1 

• 

(4 

.18a) 

^1 

= 

h + 

1^12 

A22 


2l 

=. [A22 

'i ■ 

- [A~^ 

^ 22 

A21 

] X^ 

(4 

.18b) 


with initial conditions X^ (0) = 


^ * [Agj] ai (4.19: 

with 3-0^°^ “"—I 

4.4 MODEL SIMPLIEIOATIOIT USITTG SIFG-TILAR PERTURBATION THEOTtY 
In this section, a systematic procedure for simpli- 
fying the dynamic model (4.9) of a power system using the 
theory in the preceding section is given. Since the model 
( 4 . 9 ) of any physical power system is hardly given in the 
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required form (4.10), the first and important step is to 
rewrite the given system equations in the singular perturba- 
tion form. This needs proper identification of the perturba- 
tion parameter ’e’. In physical systems ( it as well 
applies to power systems also ) ,the use of e is symbolic and 
represents the presence of a fast as well as a slow phenomena 
in the given system equations. This is primarily ascertained 
if the ratio 


Max I X. I 
i ^ 

± » i 

Min I X. I 
i ^ 

where X is the eigenvalue of [A], The next step is the 

T T T 

partitioning of y; as [ x , _z ] i.e. the slow and the fast 
varying state variables. To start with, the state variables 
associated with small time constants, small rotating masses 
and large gains are included in the group _z by expressing 
them as T^ = eT^^, = eM^ and = K^/ e, where T^, 

and are the known coefficients [46]. Then the validity 
of grouping or of the resulting reduced system (4.13a) is 
verified by the condition (4.20) given below. If the 
grouping is not valid, the state variables have to be 
regrouped and the procedure is repeated. 


Vidyasagar [60) , by considering the eigenvalues 

[ A ] as the perturbation of the eigenvalues of [Aj^] and 

[A 22 A] showed that if 

Min [Re ^^.) 

■! / I - ! 'i, -f 


of 


(4.20) 


,m-] 


2 




• » • 
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■where cr . = eigenvalue of [A„] 

J K 

= eigenvalue of [A 22 /el 

then ’m-j’ eigenvalues of [A^^] and ’ 1112 ’ eigenvalues of 

[A 22 / el'^re close to the corresponding eigenvalues of IaJ . 

We ■will use (4.20) to verify the validity of grouping the 

T IT 

state variables ^ as [x , z ] . 


4.4.1 Numerical Example to Illustrate the Grouping of 
State Variables: 

We now illustrate the above aspect of grouping the 
state variables with a small numerical example. Consider 
a single machine infinite bus example (Eig. 4.2). The 
synchronous machine model as given in Appendix 0 is consi- 
dered and the control equipments are ignored for the sake 
of brevity .Eor H = 3.5, E = 0.005 and other parameters given 
Section 4.5.1, the system equations are given by 


“a 


Aw 


A*6 


ae’ 

! 


^ It 



- 1 .58 
27.86 
0.0 
0.23 
22.82 

5.35 


0.0 

-0.06 

-0,01 

0 . 66 

-0.01 


ae; 

0.22 -148.84 

93 . 12 - 127.86 

97.12 


Aw 

1.0 

0.0 

0.0 

0.0 

0.0 


A6 

ae; 

0.0 

0.46 ■ 

-14.94 

0 

0 

2.42 


0.0 

-2.258 

- 0.17 - 24.24 

- 0.17 

j 


ae; 

0.0 

10.68 

55.52 

5.35 

-26.43 ■ 


_ae;_ 


0.0735 0.0 


in 


0.0 

0.0 

0.0 

0.0 

0.0 


44.8786 


r “1 

0.0 



0.0 

j 


0.0 

. i 

0.0 



( 4 . 21 ) 
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In (4.21), the state variables Ae" and Ae” are 

1 d 

associated with the time constants and respectively. 

Similarly AE’ and AeI are associated with the time constants 
hl ^ 

and respectively. Typical values of T^^ and t”^ for 

synchronous machines are veiy small and are in the range 

0.01 - 0.1 sec. whereas the values of f’ and T' are in 

qo do 

the range 1-6 secs, and above. Therefore we group as 


/ = [AE^, Aw, A6, AE^l and = [AE^, AS^] 


Such a separation of state variables gives the following 
matrices 1 corresponding to the slow and 

fast subsystems (4.13) and (4.14) respectively as given 
b el owl. 


[Aj^l = 


-0.9638 

-209.3901 

0.0 


1.4048 


0,0 -0.1519 -0.0310 

-0,2244 -99.2697 303.5656 

1.0 0.0 0.0 

0.0 1.3793 -9.8533 


[ ^ 22 .^ ^ 1 = 


-24.34 -0.1765 

5.3524 -26.4377 
The eigenvalues of [Al, [A^^l and are as follows: 


-2.5189+010.4483, -0.519, .-24.6926, -33.1087, 


-4.1741 


-1.7228+03.8404, -7.0425, -0.557 

-25.7830, -24.9947 


Eor the above grouping, 
satisfactory and the re 


s = 1.4 hence the grouping is 
salting reduced system..is correct. 
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T ! 

Another partitioning of ^ as x = 

®d» ^ (4.21 ) results in [ Aj^ ] and [ ^ 22 ^ 

matrices which have the following eigenvalues. 

-0.5684, -0.3109+D7.5029 

-7.3674, -24.3674, -33.5953 

for this grouping s = 0.436 and the condition (4.20) is 
not satisfied. Therefore , such a grouping will not give the 
correct reduced model. 

After grouping the state variables, the parameter e 
is assigned a small value < 1, usually equal to the smallest 
time constant and all other time constants associated with 
the decay of variables z are expressed as eT^. Equation (4 .9) 
is now in the form (4.'l0) by taking e to the left hand 
side along with z. 

4.4.2 Solution Procedure: 

(a) Reduced Model : The reduced model of (4.9) is given by 

= lAj^l + [%1 li (4.22) 

where [k^] and [B^] are defined in equation (4.13a) and the 
dimension of X^ is m, < m. Xo 

solution of X. Zq is given from the algebraic equation 
(4.1 3b). 

(b) Zeroth Order Solution of z: The zeroth order solution 
to z is obtained by solving equations (4.13b) and (4.14) 
and adding the responses as indicated in (4.15). 
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(c) Pirst Order Solution of x : The first order solution of 
X is obtained by adding the term e[X^(t) + Pq(t) 1 to 
where and x-| are obtained by solving (4.17) and (4.18a). 

(d) First Order Solution of z: The first order solution 

of z is obtained by adding the term e[Z^(t) + (x ) ] to 
the zeroth order solution of z in the step (b). (t) 

and q-i(T) are obtained by solving (4.1Sb) and (4.19) 
respectively. Higher order responses can be obtained in a 
similar manner, if necessary. 

In the above solution procedure the equations (4. 22) 
and (4.1Sa) which correspond to the X^ and X-] (slow subsystem) 
and the equations ( 4 . 14), (4.17) and (4.19) corresponding to 
q , p and q^ (fast subsystem) are solved separately and in 
different time scales. This important feature of singular 
perturbation theory allows larger time step size in the 
numerical integration because the matrices (slow 

subsystem) and t -^22 ^ (fast subsystem) are not as stiff as 
the original system matrix [A ]. This leads to computational 
advantags as shown in the numerical examples given in the 
next section. 

Thus, if an approximate response to x is to be obtain- 
ed, we use the reduced system (4.22) in step (a) only. This 
is same as the simplified model for final period in the 
state variable grouping technique t39] or the classical 
method of reduction 144] . However, if the effect of the 
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fast phenomena on the response of x is to he considered, the 
additional steps (h) and (c) are carried out. Generally, the 
first order solution of x to effect the improvement in 
response is adequate. 

4.5 NTMERim EXAMPLES 

4.5*1 Single Machine Infinite Bus System: 

The single line diagram of the system is sho'wn in 
Eig. 4.2. The system data is given below. 

Synchronous Machine : MVA = 500, H =10 sees, L = 0.03 , r^^ = 0.0, 
Xj = 1.99, x^ = 1.87, x^ = 0.3, = 0.45, x^ = = x’'= 0.21, 

’'all = °-^5, = 5.4, = 1.33, 15 , = 0.053, T"„ = 0.061. 

Exciter-Voltage Regulator System: = 50 , T ^1 - 0.1 , 

^A 2 > % = 0 . 0 , Tj; = 1 . 2 , = 0.04, T^ = 0.5, 

Tj^ = 0.05. 

Turbine-Governor System : Kg = 20, T^ =0.0,6, 

T^ = 1.3, T^ = 0 . 0 , = 0 . 5 , K 2 = 1 , = 1 , 314 rad/sec 

Transmission line 0 . 01 , xj, = 0.1 

Power transmitted = 0 . 4 - 30.3 p.u. 

The order of the full scale model is 14. The perturbation 
parameter e is identified in small time constants 
^qo» ^A 2 » % ^ 0 * ^ partitioned as 

X® = [AE’ , Am , Ad , AE', AE.^, 7 ^ 2 * ^ eV ^tV ^t 2 3 

— q U 

yslf n 

^ =[AE^q^ , A E^, 7e4» ^e 5 ’ ^e5 * 



Ill 



Transmission Infinite 
generator line bus 


Fig.A-:^ Single mac. line infinite bus systenr 
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The ratio in equation (4.20) for the above system is 
s = 1 ,1933 > 1 


Hence the partitioning is satisfactory. The initial condi- 
tions of the state variables are x(0) = 0 and z(0) = 0 . 

A step input equal to = 0.1 is applied and the 

responses are obtained for (i) the full scale system, 

(ii) the reduced system (zeroth order solution) and 

(iii) improved responses (first order solution) for both 
X and z. These responses of A 6 and AE^ belonging to the 
group X are shown in Figs. 4.3 and 4.4 respectively. The 
comments on the responses are given at the end of this 
section. 

4.5.2 Three Machine System: 

A single line diagram of a three machine system is 
given in Fig. 4.5. The line parameters and the loading condi- 
tions of the system are shown in the single line diagram. The 
machine parameters are: 


Machines 1 and 3 : 500 MVA, H = 10 secs., F _ 0.03, 0,0, 

Xj = 1.99, = 1.S7, xj = 0.3, = 0.45, = x" = 0.21, 

0.15, = 5.4, I' = 1.3, = 0.053, I" = 0.061. 


a 2. 


Machine 2; MVA = 500, H = 10 secs, 1 = 0.035, - 2, x^^-l .3, 


x^ = 0.27, x^ = 0.45, \ Xg^ = 0.15, 

T’ = 4.3, T’ = 1.4, = 0.03, = 0.04, = 0.0. 



Full scale system response 
Reduced system response 
Improved response 



Fig. 4^ Responses of state variable A6 (Single machine infinite 
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Each, ina chins is assumed to have an sxci'ter'-^regulatoi’ and 
a governor system and the parameters are same as in the 
single machine infinite bus example. 


The order of the system matrix is 42. Eive state 
variables of each maohine associated with small time constants 
as given in the single machine infinite bus example are 
grouped in z. Therefore, the dimensions of x and £ are 
27 and 15 respectively* This grouping for the given system 
gives 

s = 1 ,0762 > 1 

Hence the partitioning is satisfactory. 

The response of the system for a step input 
= 0,1 is computed for (i) the full scale system, 

(ii) the reduced system (zeroth order solution) and 

(iii) improved response (first order solution), As an 
illustration, the responses for the state variables, Aw^ 
and Ad 2 belonging to x are shown in Eigs. 4.6 and 4*7 
respectively. The response y^^ of machine 1 belonging to 
^ is shown in Eig. 4*8. All these responses are plotted 
for a period of 8 seconds. 

4.6 DISOUSSIOH 

4.6.1 Responses for Single and Three Machine Systems: 

(i) The response of the reduced system is very 
approximate in both the cases. It matches for large t, 











•whereas for the transient period of interest, there is 
significant error. 

(ii) predicted "by the theory, the first order 
solution represents significant^ improvement over the 
reduced system response. 

(iii) In I'ig.4.S,it may he noticed there is a differ- 
ence in response near t = 0 of the variable y^^ of machine 1 
(belonging to z ) . This is due to the difference in initial 
conditions of z and Z (i.e. a (0) ^ 0). This difference in 
response is known as the boundary layer jump [50]. 

4.6.2 Computational Times: 

In the numerical integration of the original system 
as well as in the singular perturbation approach, we have 
used BKG method of integration and the integration is done 
with the maximum allowable step size in each case. In order 
to compare the permissible time step sizes in the above two 
approaches, the eigenvalues of the matrices [A] (original 
system), [^j^l (slow subsystem) and [A 22 ] (fast subsystem) 
are given in liable 4.1. 

As a rule of thumb, the permissible step size in 
an explicit method is given by 

Max" lap 

i 

where is the eigenvalue of the matrix of the system to 
be integrated. Thus, the eigenvalues in Table 4.1 indicate 
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i:ha.t in "the integration of the slow subsystems involving 
equations (4»13) snd (4«18), we can use as high as five 
times and in the integration of the fast subsystems 
involving equations ( 4 . 14) and (4. 17), we can use ten times 
larger time step size than for the origial system equations. 
The computational saving in the singular perturbation 
approach is mainly derived from this fact. The permissible 
time step sizes in the integration and the computation times 
for the two numerical examples given in the preceding 
sections are summarised in Table 4.2. 


Table 4.1 : Eigenvalues. 


System 

Matrix 

Eigenvalue with 

Eigenvalue with 


max. magnitude 

min. magnitude 


[A] 

-33.408 

-0.410+30.8241 

Single machine 
infinite bus 

[Aj^] 

-8.1 

- 0 . 458 + 30.32 

-1.667 


[k22^ 

-2.57 


e =0.1 




[k] 

-46.5972 

-0.0237+30.1165 

Three machine 
system 


-5.9926+34.14 

- 0 . 02586 + 30.122 



-4.6192 

-1 .666 


e t=o . 1 
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Table 4.2 : Oomputation Times. 


System 

Integration of the 
original equations 

Singular perturbation approach 


i. llowable 
step size 

Time 

Allowable 
step size 

Time for 
zeroth 
order 
solution 

Time for 
first 
order 
solution 

Single 

0.01 

15.30 

0.05 

1.9 

6.01 

machine 

infinite 

bus 

sec. 

secs . 

sec. 

secs. 

secs. 

(inclusive of 

1 .9 secs. ) 

Three . . 

machine 

0.01 

S5 secs. 

0.05 

18 secs. 

4b secs. 

system 

sec. 


sec. 

(inclusive 
of compu- 
tinglAj^l) 

(inclusive 
of 18 secs, 
for zeroth I 
order solu- ! 
tion) 


Note: The times given in the table are GPU times on Dec System 10. 


The computation times in the table indicate that the 
zeroth order solution can be obtained in far less time as compared 
to the time required for the original system response. The first 
order solution, v^hich is a very good approximation to the response 
of the original system is obtained v/ithin 50/» of the time required 
for the original system response. VTe expect similar results vith 
other explicit methods of integration, i/^bere numerical instability 
is the dominant factor in determining the time step size. Since 
the basic motivation involved in the chapter is that of tb 
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dynamic simplification. We have not compared the times using 
the different methods of integration. 

In comparison "with the conventional modal methods, the 
singular perturbation tedinique is simple and does not need 
the computation of the eigenvectors and reciprocal "basis 
vectors of the system matrix. Although the grouping of the 
state variable is intuitive to start with, the validity of 
the reduced model is verified by the criterion (4.20). This 
needs the computation of eigenvalues of I -^ 22 ^ eigen- 

value with the largest magnitude of’ lAj^] . Once the equations 
are put in the singular perturbation form correctly, the subse- 
quent analysts is straightforward. 

4.7 COTTOLUSIOW 

In this chapter we have presented a new technique 
based on singular perturbation theory for the dynamic simpli- 
fication of the linear models of power systems. The techni- 
que is validated on two power system examples and its merits 
are brought out. The technique is simple and has potential 
for application in the dynamic studies of large scale power 
systems. 



CHAPTER 5 


L/YHAl''iIC SIMPLIPICATION USING SINGULAR PERTURBATION THEORY 
(TRiNSTENT STABILITY STUjJIES) 

5.1 INTRO LUCITON 

In this chapter, we shall extend the results of the 
previous chapter for transient stability studies, i.e. the 
application of sin^'pilar perturbation theory to the nonlinear 
model of a power system. The chapter is organised in a 
similar manner as Chapter 4 except that the system equations 
are nonlinear and the singular perturbation theory as applied 
to nonlinear initial value problem [50] is used. The problems 
associated with the simplification of nonlinear models using 
the method of asymptotic expansion are discussed. The 
computation of first order solution involves solution of 
differential equations whose coefficients are fimictions of 
the reduced system response. The technique is illustrated 
with a single machine infinite bus example. Extension to I 

i 

multimachine power system in principle is possible. However, i 
computational complexities grow enormously as compared to the 

i 

dynamic stability .problem of Chapter 4. I 

i 

5 • 2 SYSTEM MOuBL I 

The overall system model for' transient stability studi^' 
is obtained by assembling the differential e<iuat ions represent 
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the dynamics of tlie generating units and the algebraic eq. na- 
tions representing their interconnections. 

In general, the system model is in the following form. 

i = I £) > with 2 :( 0 ) = (5.1) 

0 = (J (y., s) (5.2) : 

where y, is the state vector of dimension 'm' and s is the 
vector consisting of non-state variables. In simulating the 
system for a given fault, ecquation (5.1) is numerically 
integrated and the network equations (5.2) are solved at the 
end of each integration step. The changes in the conf iguration^ 
of the transmission network due to initiation or clearance of 
a fault are appropriately reflected in equation (5.2)., 

In Section 5.4, where a numerical example is given, we 
shall give in detail the equations for a single machine 
infinite bus power system and the equations involved in the 
subsequent dynamic simplification process, i.e. computation | 
of the first order solution. 

5,3 lYlUMIC SIMPLIFICATION OP THE SYSTEM MOLBl 

The nonlinear initial value problem in the singular 
perturbation theoi^^ which forms the basis of the dynamic 
simplification of potrer system model given in (5.1) ^nd (5.2), 
is discussed in Appendix D. The basic steps in the simplrfx- 

cation procedure are as follows. 
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(i) Po rmula c ion of the given system eciuations ia the 
singular perturbation form ( equation (5.3)). 

(ii) Computa'cion of the reduced system and the zeroth 
order solution. 

(iii) Computation of the first order solution. 

In step (i) , equation (5.1) is transformed in the 
singular perturbation form as 


k - 1 ( 2 » 25 ^5 ®) 
= £ (2^5 z, e, s) 


(5.3) 


with initial conditions x(0) = 2 ( 0 ) = z° and the 

network equation (5.2) is written in the form 


G (Xj Zj s) = 0 (5.4) 

The dimensions of x and z are 21 ^ and respectively such 
that m-, + m„ = m. e is the perturbation parameter, small 

-L d, I 

and > 0 . 

I 

i 

The above process of writing the equations in the | 

T 

singular perturbation form requires, the partitioning of 2 I i 

I 

as [x^, z^3 i.e. the slow and fast state variables respec- | 

i 

tively. The physical significance of e, partitioning of the | 
state variables and the validity of the resulting reduced I 

I 

system have been uiscusssd for linear time invariant syscems i 

j 

in Section 4.3. As far as the grouping of the state variablej 
is concerned, the same procedure as discussed in Chapter 4 | 
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may be used by linearizing the given eq_uations around the 
pre-fault operating point, ll'tematively the physical basis 
for separating the fast and the slow variables may be used 
based on prior experience. In addition to tliis we have to 
verify the hypothesis (1.5) given in Appendix h, i.e. in 
equatioii (5.3) by setting e = 0, z should be expressible 
as a function of x, 

Then the resulting reduced system is given by 

lo = A s) 

= *0 (Xo. a) (5.5) 

Zo = ♦(Xo-h (5.6) 

with initial conditions X.(0) = along with the network 
equation 

i do’ S.) = Q. ^5.7) ; 

Solution of equations (5.5)-(5.7) gives (X^(t) , ^^(t)). The 
response X^Ct) is termed as the zeroth order approximation i 
to X. In physical terms, setting e = 0 to obtain the reduce(i 
system amounts to deletion of fast phenomena in the original 
model. Therefore, the response of the reduced system will, i 

i 

in general, be approximate. I 
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First Order Solution : In order to account for the effect 
of the fast phenomena on the system response, the first 
order solution of the system equations (5.3) and (5.4) is 
derived. I'he first order solution is defined as 


x(^) = ->■ £ [Ii(-fc) + 2 q(t)] + 0 (s") 

z(t) = Z^{t) g^^(T) + £ [Z^{i) + i^(T)j + 0 (e^) 


( 5 . 8 ) 


This needs the solution of the sets of differential equa- 
tions for and as given in (i).8) to (l).12) 

of Appendix n respectively. It is important to note that 
the e-quations (n.lO) and (h.ll) have coefficients dependent 

reduced system response. In 
power system problems, the analytical expression for the 
reduced system response i.e. (^^(t), ^^(t)) is rarely 
available and it lias to be stored to evaluate the elements 


of matrices in (d.lO) and (n.ll). TliLs fact makes the first 
order solution computationally -unattractive for large scale 
power systems. 

5.4 SINGLE MACIIENE INFINITE BUS EXAtlPlE 

In this section, the procedure for dynamic simplifica- 
tion described in the previous section is illustrated with | 
the help of a single machine infinite bus power system | 

example (Fig. 5.1). In the example we consider a three | 

winding model for 'bhe synchronous machine as shown in Fig. 5.| 
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Fig. 5-1 Single machine infinite bus system 
with a double circuit transmission 
line 
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and the simplified exciter-voltage regulator and turbine- 
governor models as shown in jPigs. 5*3 and 5.4 respectively. 
I'he time constants and in the above control 

circuit!, are identified as .:hown in Pigs. 5.3 and 5.4. 
Generally the time constants and are very small as 

/s. 

compared to Tg and 1^^. Therefore, by expressing them as 
and e(C the system equation are written in the required 
singular pertiurbation form (5.3). In the system equations 
the different variables used are as follows 

Xi = X2 = 0 ), X3 = 6, x^ = 

The variables and identified as shown in Pigs. 

5.3 and 5.4 respectively. Initial conditions of the state 
variables are x(0) = x'^ and z(0) = z°. 


( a) Pull scale system : The original system equations in 
the singular perturbation form are given below. 


= [-^1 + ^4 - 

^2 = - I’d - '^2 - 

where 

Pm = *5 + la ^ 2 ^ 

\ ^ - ^d) 

= X2 CO Q 

X4 = [-1% X4 + 2 i]Ae 


(5.9) 

( 5 . 10 ) 


(5.11) 


( 5 . 12 ) 
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Fig. 5*3 Simplified model of an exciter -voltage 
regulator system 



Fig. 5-4 Simtified model of a turbine- 
-governor system 
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% - [^2 + 22 ) - ^5]/(Vii) (5.13) 

£Zl = i-^x h - Vy^A (5.14) 

A 2 

£ 2 . = [-22 + K (“r - (5.15) 

^ o 0 


The network reference frame is assumed to coincide with that 
of the infinite bus. Using the transformation given in 
equation (4.4), the network equations are written as 



The terminal voltage of the machine is given by 



Simulation of the original system involves the numerical 
integration of equations (5.9) to (5.15) and the solution 
(5.16) and (5.17) at the end of each integration step. The 
fault condition is appropriately reflected in equations 
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(5. 

16) and 

(5 

.17). 






(b) 

Reduced 

System ; 

Setting e = 0 and substituting 

for 


^1 

and Z 2 

fro 

'm (5.: 

14 ) 

- ( 5 . 15 ) in equations (5 

.10), (5. 

12) 

and 

(5. 

13 ) , we 

ge 

it the 

re 

ducec system as foil ows . 

The reduced 

sys 

item variables are 

denoted as I^^ 2 ’ •* 

• and 


[Zqis Zq2] 

«> 








t)i 


[“^^1 

'1- 

\)4 - Id]/^do 


(5. 

18) 


t)2 



^e 

J 1 

3 

I 

CVJ 

■ 

1 


(5. 

19) 

where 










== 

^5 


[llmo Eg (“r - 





^e 


^1 ^q 

-1- 

I'-d ■’■(i ~ 





^3 

= 


„ 0) ) 

0 


(5. 

20) 


^4 

= 

["% 

^4 

V (V^f - vpjAj, 


(5. 

21) 


\)5 

r= 

[“^5 

-h 

h (ho + e; - ;to2» 

]/(VP 

(5. 

,22) 


The network equations for the reduced system are same 
as given in (5. IS) and (5.17) except thaf the variables x^^'s 
are replaced by 

Zq^ and Zq 2 are obtained from equations (5.14) and 
( 5 . 15 ) as 

% = <’'ref - ^t) 

(5.23) 

Z02 = K^(“r-V 

Initial conditions for are = x°, i = 1 5. 
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Structurally, the equations (5. 18) -( 5. 22) correspond to the 
set of equations (5.5) and (5.23) correspond to (5.6). 


(c) Fi rst . Order Solution : The equations for deriving the 

first order solu'.ion of = [q^^, ig obtained by 

solving the aiflerential equations co;rres ponding to (h.S) 
of Appendix au given below. 


dx 


1 

0 


■ ‘^01 

O 

X — ii 


0 

- ' 
1' 

•‘•c 


^ ‘^02 


with initial conuitions qQ-^(O) = ^\ef " V^Co"^))'*’ 

q02(0) = -Kg - Xjlob) + 

and T =: t/s. '/it!! the initiation of an electrical fault, 
the terminal voltage of the generator changes instantaneously 
ana this provicos the initial condition for (5.24). The 
differential equations for = IIPox» ••• ^05^. corresponding 
to (u.B) of Appontii:: u are as follows 


ciPQi/ci't 

, 

rj 0 

CM 

O 


0 k^^/H 

dpQ^/dt ^ 


0 0 

dPo4/dT ^ 


^ lAjs 0 

ciP05/^'f 


0 k2/(%\) 

IMWM. 




102 ^) 


(5.25) 
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with initial conaitions = 0. 

The variables = [x^ a^a = [Zu.Zij^' 

are ob^ainou rroi.i aolving the equations (5.26) and (5.27) 
given below (rei'cr to equ tions (i..l0) a;-.d (n.ll) of 
Appenuix D) . 



with initial conui tions 1 = 5, where 

j;)^(u) ia obtaint:.- iroi.i equation (5.25). 
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In the above equations, are the co:.:res ponding elements 

J- J 


of the Jacobian 


(X,. 




j and b . . ' s are the elements 
^ J 


3£o -,1 

of ], wYiere = [^'q 3 _> ^02’ ••*» ^Oml^ 


-0 " L-oi’ ^02- 

.T ... .... 


So = [«Qi* Sq 2 » •••» Som2^ obtained fi-omR.H.S. in (5.3) 

by setting e = 0. ‘Ihen the functional expressions for 


and b^j are as (aven belo-w 




'11 

II 

<»j 


s ^ ^01 

■13 

- SIT^ 



•21 

*- 3x3^ '• 


[-1 - (x^ - x^) 


■• = [-(^d - 

? r V !h T ih '“''4 ■ 

==[“In - aT - ^d 3^ 


(5.28) 


(5.29) 


- ^1 9X^ -^d 3 x^ 


■Cl § X-, 


(5.30) 
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3f, 


'23 


fcn T — 


02 




- r_x iia . I tia 


d 3 X. 


■'ll 


'13 


9^1 

3X^ 

3x^ 


h 

nC' 3 X-| 


K, 3 ?, 

A t 


3X: 




3 Id (V^a' 


■q 3 X, 


]/M 

(5.31) 

(5.32) 


(5.33) 


In the above expressions, the partial derivatives of Iq^ 
anu V,j. w.r.t. x^ and x^ are computed from the network equa- 
tions (5.16) and (5.17) and are dependent on X^(t) . Hence 

a. -’s and b. .'s are, in general, ftinctions of (^^(■t) s ^^(t ) ) . 
3-3 3. 3 

1' 

Equation (5.26) is solved by eliminating 
by using (5.27). is a known response derived from 

equation (5.13) as given below. 

dY, 

2oi = dt“ 

(5.34) 


Zo 2 = -^g ^2 

= -s' LV^e"^ ^^2“ 

§ 

In (5.34) -we have 


lit . ^ (5.35) 

•gfr" 3X^' dt 3X^, dt 

which is evaluated along ^^(t). Wirh odXLS.Sbep we can 
integrate C5.26) to oMato Xi(t) roi hence Z^{t) from (5.27) 
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To complete the first order solution, it now remains 
to solve for whose differential equations 

corresponding to (D,l2) of Appendix D are given by 


dill" 

dT 


1 

1 

> 

0 


^11 

^^12 
d t 


0 

_ i_ 

<A. 


‘I12 

— 


with Initial conditions 

<lll(0) = -Z^J^(O) 

■312(0) = -Zi2{0) 

In the present example, the equations for the 
boundary layer terms q^, and £2 given in (5«24)> 

(5.25) and (5»56) respectively are linear and time 
invariant. In general, these equations have the form 
as given in (3). 8), (D.9) and (D.12) of Appendix 1. It 
is important to note that in the boundary layer equations 
the variables (Z^, Z^) are assumed to be constant at their 
initial values (lo(0), Zq( 0)), This assumption is justified 
on the basis that 3.^, 2o % represent the fast variables 
and decay rapidly. ^^(t)), being slow variables, can 

be assumed to be constant at their initial values (^q( 0 ) s'^q( 0) ) 
in the region near t = 0 (boundary layer region). 

The first order solution of (x, z) is obtained by 
adding the responses as given below. 
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To complete the first order solution, it now remains 
to solve for = [q -,-, , whose differential eqnati 

corresponding to (D«l2) of Appendix D are given hy 


ons 


■ daii] 

r* 

dT 


dqi2 


_ dT 




0 




(5.36) 


with initial conditions 


qil(O) - -2^1(0) 

In the present example, the eqmtions for the 
"boundary layer terms 3.^, 2© 3*1 1^ (5.24)» 

(5.25) and (5.36) respectively are linear and time 
invariant. In general, these equations have the form 
as given in (D.8), (D.9) and (D.12) of Appendix D. It 
is important to note that in the "boundary layer equations 
the variables (^q» ^ 0^ are assumed to be constant at their 
initial values (Xq( 0), 2^(0)). This assumption is justified 
on the basis that 3.^, 2^ and represent the fast variables 
and decay rapidly. (Io(t), 2 Q(t)), being slow variables, can 
be assumed to be constant at their initial values (^q(O) ,^q( 0)) 
in the region near t = 0 (boundary layer region). 

The first order solution of (x, z) is obtained by 
adding the responses as given below. 
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Hesponaea Jor the above disturbance are plotted for 
(i) the full 6ca.'.e system, (ii) the reauced system and (iii) 
the improved response (first order solution). As an illustra- 
tion, the responses for the state variables 6, <o and are 
■shown in Pics. 5.5, 5-6 and 5.7 respectively. 

5.5. A-ix. COHCLUSIOW 

Ihe reopontjcfj in Pigs. 5. 5-5.7 clearly indicate that the 
deletion of two ar.-ai time constants and 1'^ results in the 
rt'uuccu nystem response differing fI^Dm that of the original 
£<yr. toi.'i rooponsc. me effect of these tiiae constants on the 
ays tow responoe ir; retrieved to a sufficient accuracy by 
uorivin^^, the first order solution. It is seen that first order 
tippioximat ion improves the response of the reduced system 
conaiuerably. However, this improvement in response is at the 

cost of increased computation time. In a large scale power 
uyateni model, deletion of small time constants as in the clas- 
sical method [44j I’liay lead to erroneous response. The singilLar 
perturbation tGclmiq.ue proposed in this cha^pter can be used to 
obtain more accm‘ate system, response. However, as is evident 
from the numerical example, the major problems in obtaining 
improved responses are: (i) separation of the slow and fast 
variables, i.e. writing the given system equations in the 
required form (5.5) <3nd (ii) solution of the differential 
equations with vai' table coefficients, Ihe computational burden 
associated with tl'.e technique increases rapidly with the size 
of the system. Hence before it can become a practical tool for 
inaustry, improved algorithms for quick solution of the differ- 
ential equations involved (with variable coefficients) are needed. 
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CHAPm 6 


COICLUSION 


6.1 summary All^ CONCLUSIOK 

In this chapter, we summarise the important results 
of Uiis thesis and indicate suggestions for future work. 

(j^) Identif ica tion of Coherency : Among the existing 
mctl'.oas for identifying coherency, the linear simulation 
method proposed by Podmore [24] is the most promising one, 
But it involves nuraerical integration and its repetition 
for each fault location in the system. On the other hand 
it ia convenient for a power engineer to deal with an 
inoex representing the coherent behaviour between two 
machines rather than computing and processing the swing 


curves . 

With this background, we have proposed in this thesis 
two methods for identifying coherency. 

(i) A u.sthod involving the comparison of swing curves 
obtained from the closed form solution. 

(ii) fhe direct method using the concept of coherency 
index reflecting the coherent behaviour between two 


machines. 
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ihe analyt/ical approach losed in the above methods, 
eliminates the numerical integration of the system eq.uations 
anu provides a better understanding of the system dynamics 
from the coherency viewpoint. Repetitive computations for 
.nul ciple fault locations are also eliminated by the appro- 
priate use of linear transformation. Therefore, the tech- 
nique is suitable for identif3ring (i) the coherent groups 
for a set of faiiLt locations in the given study system, or 
(ii) the coherent groups for different study systems. The 
proposed methods have been validated on realistic power 
systems. For large scale power systems, computation of all 
eigenvectors anu reciprocal basis vectors may involve 
increased computer time. However, the use of the concept of 
dominant modes may alleviate tliis difficulty, 

( b) i^ecompositi o n of Power Systems and Construction of 
Coh e rency Based Equivalents : It is logical to represent 

the dynamics of those generators in detail, which are 
.lOnsitive to the given fault than those Wiiich are insensi- 
tive to it. The prevalent basis for modelling of generators 
does not consider tliis aspect of electromechanical inter- 
action between the machines. The decomposition technique 
proposed in this thesis incorporates this feature and 
divluGS the power system into different regions for the 
purpose of modelling the generators. Since, decomposi'cion 
ana coherency analysis are interrelated, a procedure has 
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been proposed to ensure that a coherent group is not spread 
over tjc regions. Such a decomposition facilitates the 
; iouping of colierent generators in each region separately 
nnei hcr^o simplifies the grouping procedure by reducing the 
cciMjiurioons . For rhe sake of completeness of the topiCj 
r: VO also presented the case study for a realistic power 
tijii uua involving, ue composition, coherency analysis and 
i-he nubnequen i- ag, .rogation of coherent generators. 

i c i jimplific ation Using the Singular Perturbation 

i'hoor.v ( uynoraic Stability Studies) ; Ihe conventional modal 
moi-houa for simplifying linear models of power systems 
require cociputution of the eigenvectors and reciprocal basis 
Yoclora of the system matrix, -which may be very large in 
iiizo. The state variable grouping technique [39] and the 
cluaiiical method [44] are simple but need a physical insight 
into the uynamics of ’the system to neglect certain variables 
or small time constants etc. Moreover, these methods give 
r‘ loss accurate reduced order model of the system. 

In this thesis, we have proposed a new technique 
bri.'ieu on 3ln(5ulai- perturba-tlon tlieoiy for dynamlo simplifi- 
cniioti of linoar power system models. It involves separa- 
tion of the linear dynamics of the given system into slow 
anu rnoi subsystous. which are integrated separately in two 

These responses are added to obtain 


ui If ©rent time ocales. 
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the 2’esponse of die original 4 .. 

-i-^-max system. Separation of the 

vector into slov and fast variables is achieved 
U’.rouc.h a systematic procednre, vMoh ensures the validity 
or the resultlnr, reuuoed system. Ihe existing methods in 
the literature nedeot the fast phenomena in the system 
.oltof^itlier. In the proposed teohnlaue, this is taken into 
coiu;iUfra(.ion to obtain an improved response through first 
orucr nolution. fbe technique gives computational advantages 
o\ur ihe integiViwion of the original system equations. Ihe 
wothod haa been valioated on two power system examples, 
uosuputiiiion of 'iiic first order solution requires higher 
computer ntoraj';^© and therefore, it has to be supported by 
efficient prograffii.iing methods. 


^ wynaiuic Sim pl ifi c ation Using Singular Perturbation 
f oclmique (franaient stability Problem) ; Singular pertur- 
bation theory has been extended to simplify the dynamic 
raodoln for transient stability studies. As compared to the 
classical method of reduction, where small time constants, 
rotating masses etc, are totally neglected, the proposed 
mothod takes their effect into consideration in the first 
order solution. Hov;ever, as compared to the linear problem 
of Chapter 4, it entails the solution of the , differential 
equations with variable coefficients, 'fne technique has 
boen demonstrated using a single machine infinite bus example. 
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ii3Ci<enGion of "the 'i/ecli3iic[ue to multi-inacliine system, in 
principle, is possible but needs further refinement to make 
11 computationally attractive. 


G.2 SCOPE B’OR I'Ullf HER RESEARCH 


(i) In the coherency analysis, computations involving 
eigenvectors and reciprocal basis vectors could be reduced 
ii Vi'e conuiaer only the dominant modes in the closed form 
solution. Our experience with different systems indicate 
tluit it is sufficient to consider the modes in the lower 
natural frequency range for identifying coherency. This may 
be uue to the fact that the modes of lower frequency are 
dominant [76]. I’Ms assumption need to be investigated in 
uotnil. Then the corresponding eigenvectors and reciprocal 
bniiiti vectors coulci be easily derived by using the modified 
l)ov.'or methou given in Eef. [77]. This will make the direct 
methoa of coherency analysis more efficient. 


(ii) With the assumption of linearity in coherency 
an.nlysi=, the generators near a fault location are sometimes 
v,ro.w-;ly ereuped as coherent. Therefore, Instead of assuming 
a linear model for the complete system, it is logical to 

represent the study system by a nonlinear model and the 

n.;v,raov. mndftl The coherency analysis 
external system by a linear model. 

■ T A atroi nniiient mav lead to useful results, 
with such mathema'ijical devel p 
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(iii) The existing techniques for dynamic aggregation 
of coiiierent generating units use frequency response approach. 
The time domain approaches could also he considered or the 
problon; be analysed as the parameter identification problem 
of moaern control theory. 

(iv) i/ynamic simplification of multi -machine system, 
for ~ rjinfJient stability studies using singular perturbation 
theory results in enormous computational complexities and 
need further rcrincment. 

(v) In Ref. [57], singular perturbation approach has 
been used to retrieve the high frequency oscillations occur- 
ring between the coherent machines. For large scale systems 
thin needs to be examined further using criteria for coherency 
uo well singular perturbation approach. 
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appbtdix a 


1. UO'FD FORM SOLUTIOl^ OP EQUATIOI^ (2.5) 

The solution of (2.5) is given ty [ 78 ] 

/ V 

~pf ® ® 

This is expressed as the linear combination of decoupled 
modes ns 


X,t 2n-2 X.t 

Ul ^ ^ ^ i= 2 m +1 ^ 


where X^, arid are defined in equation (2,6). The first 
’ 2 ni' terms corresponds to *m’ complex pairs of eigenvalues and 
the rest are corresponding to real eigenvalues. 


The terras corresponding to a complex pair X^ and 
in (A, 2) are 

y<Y X>e'' +X -i<T (A. 3) 

=“ ,Xg> e + >_e 


Sin X 


m+1 * Ht+I 


^ and are the complex oonjy-gates of 


%+1 


and respectively, equation (A. 3) reduces to 


X^ t 


(A. 4) 


= 2 Re CX^ ® ^ 

By rearranging the terms, equation (A.4) is simplified to 

I ( OosB-it + sJl Sin git)X>i 

_ Sin 01 - ^ -1 ^ 


*SE ^ 


(A. 5) 
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The- response oontrliuted by - an' oscillatory modes Is then 
generalised from (A.5) and equation (A.2) is written as 

. I _ n ,n 2ri-2 a . t 

~ Sin 3a t - i. los 3.-t)X. + I x! ii’. e ^ 

^ i=2m+1 ^ 

(A. 6) 

whert 

V » ^ Ii» 

The rciaponae (A. 6) is cast in the matrix form by rearranging 
the terms as 

~pf “ (A. 7) 

where the columns of E S] and the elements of ^.re defined 

in equation (2*8). 


2. 7U)3ED FORM SOLUTIORT OP EQUATIOI (2.4) 

The solution of equation (2.4) is given by [781 


X, = [Blu dr 


(A. 8) 


This is expressed as a linear combination of modes as 
9m t ^^(b-T) 2n-2 t ' ^ 

if iil J ^ ^ 0 


where u = iBlu. 


Equation (A. 9) Is simplified to 
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m 

ii-r ” 'I 


t=1 


(fi + (f^k;^ - g.k;_)x^ 

2n-2 , k! a.t 

* t 4 # ^ > 

io2m+1 


vjvjf 

1 - 6 


gl - li 


a . Cl. "ti 't ^ , 

u ^(l-e ^ Cos gj^t) + 

a.t tti ctj "t 

^(1-e ^ Cos e^t) +-~(e Sin 

*^1 i 


i 


61 


K. «<y!,u>. ICj=<Ii. i>. ^ gj 


-i '1 


Ejuation (A-'O) is rearranged as 


X, _ 
~,f 


T PK^-e ^ Cos 3j_t) + (® ' Sin 2j_t) 

^ **»»X 

i ■ 1 


2n-1 


a. t 

. , £i(-1 + ® ’ 

i=2ra+1 


are defined in equation (2.6) 


where 2^ £i 

equation (A.11) is oast in the matrix form as 


= fPl 


where the clcmen-ts of Zp 
defined in (2. 6), 


(t) and the columns of iPl 


(A. 10) 


( 1 . 11 ) 

Then 

(A. 12) 

are 



appendix b 


DA!EA BOR UPSEB SYSTEM 


rnirr.bur of buses = 71 Number of machines = 13 

I'ujr.bcr of lines = 94 Number of shunt loads = 23 

Baste KVA - 200 

BUS MIA 

name Generation Load Power Btis voltage 

Uo. 

Mag. Ang. 


1 

OBliA^XTll 

2 

0B{U£:XT400 

3 

0BRA{T)11 

4 

0BRA{T)220 

5 

OBRA(T)132 

0 

0BRA{H)11 

7 

OBRA(H)132 

B 

HIHll 

9 

RIHI32 

10 

ROBI32 

11 

MUGI32 

12 

MUG '’20 

13 

MAUI32 

14 

MAU220 

15 

OORll 

16 

GCHI52 

17 

G0R220 

18 

KiiAL132 

19 

BASTI32 


630.667 

169.722 

0.0 

0.0 

506.0 

149.492 

0.0 

0.0 

0.0 

0.0 

98.0 

31.986 

0.0 

0.0 

297.0 

124.227 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0,0 

0.0 

0.0 

0.0 

184.08 

114.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

12,8 

8.3 

0.0 

0.0 

185.0 

130.0 

80.0 

50.0 

155.0 

96.0 

0.0 

0.0 

100.0 

6 2.0 

0.6 

0.0 

0.0 

0.0 

73.0 

45.5 

36.0 

22.4 

16.0 

9.9 

32.0 

19,8 


1.030 0.0 

1.057 -5.931 

1.025 -1.406 

1.054 -6.061 

1.047 -5.064 

1.025 0.397 

1.045 --5.723 

1.025 -0.204 

1.0433 -4.976 

1.026 -7.413 

1.023 -11.762 

1.003 -10.057 
1.007 -13.859 

1.004 -11.984 

1.005 -9.413 

1.024 -16.682 

1.005 -14.460 
1.009 -17. 881 
0,981 -19.786 


continued. . 
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Bus 
No . 

Name 

Gener€ 

20 

GONI32 

0.0 

21 

FAII32 

0 

0 

22 

SUL220 

0.0 

23 

SULI32 

0.0 

24 

SUL400 

0.0 

25 

;.LL132 

0.0 

26 

A1L220 

0.0 

27 

fUTll 

304.0 

20 

m;.t400 

0.0 

29 

Hi 

261.0 

30 

P/.NK132 

0,0 

31 

P/.KK220 

0.0 

32 

P;.NK400 

0.0 

33 

LUCK220 

0.0 

34 

1UCKI32 

0.0 

35 

1UCK400 

0.0 

36 

SITI32 

0.0 

37 

Sm.132 

0.0 

3B 

BARI32 

0.0 

39 

Kii/iTll 

25.0 

40 

KHATI32 

0.0 

41 

MORI32 

0.0 

42 

KOR. ’0 

0.0 

43 

M0R400 

0.0 

44 

liiuMll 

180.0 

45 

it.M132 

0.0 

46 

NC.HI32 

0.0 

47 

Mi.I220 

0,0 

48 

i Hi’.RJll 

341.0 

49 

1 H.'.RJ220 

0.0 


Generation load Power 


BiJS Voltage 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


27.0 

52.0 

0.0 

75.0 

0.0 

133.0 

0.0 


76.287 0.0 

0.0 30.0 
70.506 0.0 

0.0 120.0 
0.0 160.0 
0.0 0.0 

0.0 0.0 

0.0 112.0 
0.0 0.0 
0.0 50.0 
0.0 147.0 
0.0 93.5 

30.583 0.0 
0.0 0.0 
0.0 255.0 

C.O 0.0 
0.0 0.0 
55.037 0.0 
0.0 0.0 
0.0 78.0 


0.0 

256-0 

0.0 


234.0 
0.0 

295.0 



Mag. 

Ang. 

16.8 

1.003 

-.21.858 

19.8 

1.016 

-21. 376 

0.0 

1.008 

-17-625 

46.6 

1.036 

-19-815 

0.0 

0.985 

-16-193 

82.5 

1.039 

-17. 441 

0.0 

1.016 

-14-382 

0,0 

1.025 

- 5.455 

20.0 

1.055 

-12.342 

0.0 

1.025 

-15.519 

74.5 

1.057 

-19.947 

99.4 

1.014 

-22.817 

0.0 

1.025 

-15.936 

0.0 

0.998 

-25.515 

69.5 

1. 5O6 

-29.639 

0.0 

0.966 

-a. 803 

32.0 

0.994 

-34.030 

9 2.0 

0.967 

-57-982 

88.0 

0.972 

-57-638 

0.0 

1.025 

-34-166 

0.0 

1,018 

-36.913 

123.0 , 

1.007 

-34.778 

0.0 

0.979 

-31.484 

0.0 

0.946 

-27.835 

0.0 

1.025 

-24.448 

0.0 

1.048 

-30.514 

38.6 

I.O2O 

-32.678 

145.0 

0.989 

-29. 053 

0.0 

1.005 

-25.781 

183.0 

0.997 

-30.500 


continued. • 


163 


Bus Name Generation load Power Bus Voltage 

Ko, Mag, Ang. 


0.0 0.0 

0.0 0.0 


50 MUR220 

51 MUH132 

52 MlTRnOO 

53 r^iiR220 

54 M2.iiRl32 

55 SiUM220 

56 Si.il»220 

57 3«.11.132 

58 lUJiD132 

59 H00132 
GO r<00220 

61 I<1S220 

62 HIS132 

63 Dr; HI 32 

64 YAMIIll 

65 

66 
67 

ijB MGNll 

60 I u. II 2 20 

70 II/.LD132 

71 SH.4220 


0.0 0.0 

0,0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

0.0 0.0 

300.0 72.882 


89.0 26.689 

0.0 0.0 

0.0 0.0 

0.0 0.0 


40.0 

24.6 

227.0 

142.0 

0.0 

0.0 

0.0 

0.0 

108.0 

68.0 

25.5 

48.0 

0.0 

0.0 

55.6 

35.6 

42.0 

27.0 

57.0 

27.4 . 

0.0 

0.0 

0.0 

0.0 

40.0 

27.0 

33.2 

20.6 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

14.0 

6.5 

0.0 

0.0 

0.0 

0.0 

11.4 

7.0 

0.0 

0.0 


0.977 -32.056 
1.005 -35.509 

0.957 -26.552 

0.972 -32.845 
1.004 -35.3I8 

0.985 -31.313 

1.013 -29.253 
1.016 -29.967 

1.018 -30. 261 
1.013 -30.715 

1.008 -29.974 

1.019 -27.874 

1.044 -28.940 

1.042 -27.695 

1.025 -19.401 
1.057 -24.662 

1.025 -20.610 
1.055 -25.641 

1.025 -19.829 

1.050 -25.363 
0,998 -34.790 

0.927 -35.329 


Yi.iai220 0.0 0.0 

96.0 25.656 

y..MI,IV32 0.0 0.0 


continued. . . 


GxJlfcRATOR MTA 


r 0 n # A ‘0 

. Bua Prom 

Name 

0 

Inertia Tirans. 
oast, (sec.) React. 

Damp 

const. 

1 

1 

OBi?i.dXTll 12.900 

0.0488 

0.0 

2 

3 

CBR/ 4 T )11 

9.000 

0.0690 

0.0 

3 

6 

0 BEi.(n)ll 

1.923 

0.5278 

0.0 


6 

RIHll 

6.648 

0.2010 

0.0 

5 

15 

GORll 

2.590 

0.2440 

0.0 

0 

21 

]yL:Tll 

2.550 

0.1455 

0.0 

7 

29 

Pi.mi 

2.700 

0.1355 

0.0 

A 

39 

KHilTll 

1.035 

1.0820 

0.0 

0 

4 4 

fL'JIll 

5.430 

0.2270 

0.0 

10 

48 

E'.RJll 

7. 808 

0.0595 

0.0 

11 

64 

YAMIIll 

7.140 

0.1300 

0.0 

1? 

66 

Y 1 MI,IV 11 

2.176 

0.4860 

0.0 

13 

68 

M/.N11 

1.473 

0.5990 

0.0 



II HB D/.TA 



Linu 

i;o. 

Prom To 

Bus Bus 

line Impedance 

charge 

Turns 

Eatio 




r\ ACT rrH 

n 0 

1.05 


1 

2 

5 

4 

5 

6 

7 

8 

9 

10 

11 


9 

9 

9 

9 

10 

7 

12 

11 

14 

13 

17 


8 

7 

5 

10 

11 

10 

11 

13 

13 

16 

15 


0.0320 

0.0660 

0.0520 

0.0660 

0.0270 

0.0 

0.0600 

0.0 

0.0970 

0.0 


0.0780 

0.1600 

0.1270 

0.1610 

O.O7OO 

0.0530 
0.1430 
0.0800 
0. 2380 
O.O92O 


0.0090 

0.0047 

O.GI4O 

0.0180 

O.OO7O 

0.0 

0.0360 

0.0 

0.270 

0.0 


1.00 

1.00 

1.00 

1.00 

1.00 

0.95 

1.00 

1.00 

1.00 

1.05 


continued 
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Line 

Ho. 

From 

Bus 

To 

Bus 

12 

7 

6 

15 

7 

4 

14 

4 

3 

15 


5 

16 

4 

12 

17 

12 

14 

18 

17 

16 

19 

2 

4 

2C 

# 

26 

21 

2 

1 

22 

31 

26 

25 

26 

25 

24 

25 

23 

29 

22 

23 

26 

24 

22 

27 

22 

17 

2B 

2 

24 

29 

23 

21 

30 

21 

20 

51 

20 

19 

52 

19 

18 

53 

18 

16 

34 

28 

27 

55 

30 

29 

56 

32 

31 

37 

51 

50 

38 

28 

32 

39 

31 

33 

40 

51 

47 

41 

2 

52 

42 

33 

34 

43 

35 

33 


Line Impedance ‘^‘^charge 


Turns 

Eatio 


0.0 

0.0 

0.0 

0.0 

0,0160 

0.0160 

0.0 

0.0 

0.0190 

0.0 

0,0540 
0.0 
0. 2040 
0.0 
0.0 

0.0480 

0.0100 

0.0566 

0.0720 

0.1460 

0.0590 

0.0500 

0.0 

0.0 

0.0 

0.0 

0.005’ 

0.0150 

0.0100 

0.0158 

0.0 

0.0 


0.2220 

0.0800 

0.0530 

0.1600 

0.0790 

0.0790 

0.0800 

0.0620 

0.0950 

0.0540 

0.1670 
0.0800 
0.5200 
0.0800 
0.0840 
0. 2500 
0. 0200 
0.1412 
0.1860 
0.3740 
0.1500 
0.0755 
0.0810 
0.0610 
0.0930 
0.0800 
0.0510 

0. 0640 
0.0790 
0.1570 
0.0800 
0.0840 


0.0 

1.00 

0.0 

1.05 

0.0 

1.00 

0.0710 

1.00 

0.0710 

1.00 

0.0 

0.95 

0.0 

1.00 

0.1930 

1.00 

0.0 

1.05 

0.1500 

1.00, 

0.0 

0.95 

0.0130 

1.00 

0.0 

0.95 

0.0 

0.95 

0.0505 

1.00 

0.3353 

1.00 

0.0140 

1.00 

0.0050 

1.00 

0.0100 

1.00 

0.0C40 

1.00 

0.0080 

1.00 

0.0 

1.05 

0.0 

1.05 

0.0 

0.95 

0.0 

0.95 

0.6706 

. 1.00 

0.0580 1-00 

O.I77O 

0.5100 1‘00 

0.0 

0.95 

0.0 

0.95 


contin-ued . 
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t 


Line 

No. 

From 

Bus 

To 

Bus 

44 

35 

24 

43 

34 

36 

46 

36 

37 

47 

37 

38 

8 

40 

39 

4 0 

4 0 

38 

50 

38 

41 

51 

41 

51 

32 

4 2 

41 

53 

43 

4 2 

54 

47 

49 

55 

49 

49 

56 

49 

50 

57 

49 

42 

58 

50 

51 

59 

52 

50 

60 

50 

55 


61 50 53 

62 53 54 

63 51 54 


64 

55 

56 

65 

56 

57 

66 

5” 

59 

67 

59 

58 

68 

60 

59 

69 

53 

60 

70 

45 

44 

71 

45 

46 

72 

46 

41 

73 

46 

59 

74 

60 

61 


Line Impedance 


0,0062 

0.0790 

0.1690 

0.0840 

0,0 

0.0890 

0.1090 

0.2350 

0.0 

0.0 

0.0210 

0.0 

0.0170 

0.0370 

0.0 

0.0 

0.0290 

0.0100 

0.0 

0,0220 

0.0160 

0.0 

0.0280 

0.0480 

0.0 

0.0360 

0.0 

0.0370 

0.0830 

0.1070 

0.0160 


0.0612 
0. 2010 
0.4310 
0.1880 
O.38OO 

0. 2170 

0.1960 

0.6000 

0.0530 

0.0840 

0.1030 

0.0460 

0.0840 

0.1950 

0.0530 

0.0840 

0.1520 

0.0520 

0.0800 

0.0540 

0.0850 

0.0800 

O.O72O 

0.1240 

0.0800 

0,1840 

O.I2OO 

0.0900 

0.1540 

O.I97O 

0.0830 


0. 2012 

1.00 

0.0220 

1.00 

0. 0110 

1,00 

0.0210 

1.00 

0,0 

1.05 

0.0250 

1.00 

0.0220 

1.00 

0.0160 

1.00 

0.0 

0.95 

0.0 

0,95 

O.O92O 

1.00 

0.0 

1.05 

0.0760 

1.00 

0.0390 

1. 00 

0.0 

0.95 

0.0 

0.95 

0.0300 

1.00 

0,0390 

1.00 

0.0 

0,95 

0.0060 

1.00 

O.OI7O 

1.00 

0.0 

1.00 

0.0070 

1.00 

0.0120 

1.00 

0.0 

1.00 

0.0370 

1,00 

0.0 

1.05 

0.0100 1.00 


O.OI 7 O 1**^^ 
0.0210 1*^^ 
0.0160 l»O0 

continnBd . * 



Line 

n 0. 


Prom 

bus 


To 

Bus 


Line Impedance Turns 

cnarge patio 


75 


61 


62 


0.0 


76 

9t 

62 

0.0420 

77 

62 

63 

0.0350 

7B 

69 

68 

0.0 

79 

69 

61 

0.0230 

RO 

67 

66 

0.0 

B1 

65 

64 

0.0 

B2 

65 

56 

0.0280 

B3 

65 

61 

0.0230 

m 

65 

67 

0.0240 

a5 

67 

63 

0.0390 

B6 

61 

42 

O.O23O 

S7 

57 

67 

0.0550 

68 

45 

70 

0.1840 

89 

70 

38 

0.1650 

9C 

33 

71 

0.0570 

91 

71 

37 

0.0 

92 

45 

41 

0.1530 

93 

35 

43 

0.0131 

94 

32 

52 

0.0164 


0.0800 

0.0 

0.95 

0.1080 

0.0020 

1.00 

0.0890 

0.0090 

1.00 

0.2220 

0.0 

1.05 

0.1160 

0.1040 

1.00 

0.1880 

0.0 

1.05 

0.06 30 

0,0 

1.05 

0.1440 

0.0290 

1.00 

0.1140 

0.0240 

1.00 

0.0600 

0.0950 

1.00 

0.0990 

0.0100 

1.00 

0.2293 

0.0695 

1.00 

0. 2910 

O.OO7O 

1.00 

0.4680 

O.CI2O 

1.00 

0.4220 

0.0110 

1.00 

0.2960 

0.0590 

1.00 

0.0800 

0.0 

0.95 

O.388O 

0.0100 

1.00 

O.I3O6 

0.4293 

1.00 

0.1632 

0.5360 

1.00 



SHUHT LOAD DATA 


S.Ko. Bus No, Shunt load Admittance 


1 

2 

2 

13 

3 

20 

4 

24 

5 

28 

6 

31 

7 

32 

B 

34 

9 

35 

10 

36 

11 

37 

12 

38 

13 

41 

14 

43 

15 

46 

16 

47 

17 

50 

18 

51 

19 

52 

20 

54 

21 

57 

22 

59 

23 

21 


0.0 

-0.4275 

0.0 

0.1500 

0.0 

0.0800 

0.0 

-0.2700 

0.0 

-0.3575 

0.0 

0. 2000 

0.0 

-0.8700 

0,0 

0.2250 

0.0 

-0.3220 

0.0 

0.100 

0.0 

0.3500 

0.0 

0. 2000 

0.0 

0. 200 

0.0 

-0. 2170 

0.0 

0.1000 

0.0 

0. 3000 

0.0 

0,1000 

0.0 

0.1750 

0,0 

-0. 2700 

0.0 

0.1500 

0^ 0 

0.1000 

0.0 

0.0750 

0.0 

0.0500 



APPEMDIX G 


G,1 SYNCHRONOUS MACHINE 

Cnnsider a six winding model of a synchronous machines 
riti Shown in Pig. The flux current relationship in the 

hybrid characterisation with the assumption of common flux 
1 inking all the windings in each axis is given as [64] 

Si root Axis Characterisation ; 

If 

»> ^md 

I! ” 

1 ^-4 ) 

fk. 

^kd SL ^kd Jt ^f s- 

( 0 . 1 ) 

1 » J- + 'I 

.more ^ == x^^ ^ ^ 

%d (0.2) 

and xj = + ^a£. 
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*fl 




Fig. C-1 Multiport representation of a 
synchronojs machine 
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172 



H tchlnr 

A»airatnR that xj = the synoteonous maohlne moflel 

vilth n-spo-.t to its own fl-q axes frame and Inoluslye of 


rotor riyn'^nilcs is given T^y 
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Thf rlifferent el«^ents of the matrices are defined as 


fnllowa; 




» 


*'do(-a 
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Thf tin.-' conat-^nts and and are defined in 

th'- uainl rarinner as given in Ref. [62], 

C,? KX^ITER-VOLTAGB R15GUIAT0R STSTM 

The state space representation of the exciter-voltage 
r'g'ilatnr nystcni shown in Hg, 0.2 is given hy 


t,. « 

wh r re 

.« y^,2, 

^ ^rcf 


+ [0, 


me 


' 2 b Ve 



T 


(0.8) 


G.3 TURBIW-GOVERNOR SISTM 

The state space representation of the turbrno governor 
nyaten nhown In Fig. 0.3, which can represent both the systems 
rcr hydro and steam turbine by properly selecting the parameters, 

in given by 

it = t to,.! + !0tm>2m + it "t 

where 

yt2» ^t5^ 


u. 


A* 


‘t " ” tno 

^ ^ (r p.) nnd (C.9) can be easily derived 

The details of models (0.3) an ^ 

« n ? nnd 0 3 respectively along with tTae machine 

from the Pigs* 0.2 and re p 

. n.P not given here for the sake of brevity, 
eijuations and are no s 




Fig. C-3 Simplified linear model of a turbine- 
-governor system 












APPENDIX h 


NQH^LIlfEAR INITIAL VAIiUE PROBLEM 

i'ho non-lincar initial value problem given below is 
atiapted from Ref. [50], 

Conaiaer a nonlinear system 

k “A s) 

(D.l) 

‘-1=5 (2£» z, e) 

with x(0) = X® and z(0) = z°. The dimensions of x and z 
are and iHg respectively, 'e' is the perturbation 
pararnoter small and >0. f and £ are continuously differ- 
ontiublc functions and have a^ymptomic expansions 

f (x, z, c) V I f. (x, z) 

(D.2) 

jg {x, z, e) -V. I Ba 2 ) 
j=0 

as e >0. The reduced problem is 

k ^ £ (x, z) with x{0) = x° (I).3a) 

= Bo 


0 


(I).3b) 
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A 80 u me that 

(i) There i3 a continuously differentialole fmiction_$ ( 2 ) 
such that 

£0 = Q (U.5) 

and the resultin/!; reduced system has a (unique) solution 

'•'he sake of brevity we will not use 
the arf^ment (t) v/ith ^ and in the subsequent analysis. 

8 /^r i 

(ii) The matrix [-y^] is nonsingular and stable along 

(iii) The matrix [|~] is stable at (Xq( 0), X) for all 
value a of K between ^(0) and z°. 

The asymptotic solution of (h,l) is an additive 
i'unction of the variable 't* and a stretched variable 
1 : ~ t /£ , which tends to <» as e -»• 0. The solution of (D.l) 
io written in the form 

x(t,e) = X(t,£) + £ :2 (i:,£) 

( 1 . 6 ) 

z{t,€) = Z(t,E) + a {^,e) 

where x, X, z and Z all have asymptotic expansions given 
in equation (4,12). The outer expansion (X, Z) and the 
boundary layer expansion (e£, a) are explained in 

Section 4.3. 
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’.‘e shall obtain the zeroth and the first order solu- 
tion an defined in (4.15) and (4.16) respectively. For 
l];in purpose, the different terms in (u.6) are obtained 
i y UMin/:r laylor oories expansions of f and around 
(l^. 4^^, 0) and substituting the expansions for (I, Z) and 
i“) (.^*1), On equating the coefficients of we 

hnvu the followini, equations. 


Ho a’icou nyotern 


4 = 4 (4- V 


(D.7) 


4 


♦( 4 ) 


2^ and jg^ are given by 

“4 a (t) (D. 


8) 


with initial condition 1^(0) 


- Zo(0) 


^lo r!4 4(°^ (0.9) 

-it “ ^3 2 


with initial condition 

j^o (oo) = 0. 

Ihe terms 4 and 4 are given by 
3f 

4 = [^r ^ -1 '•32 


,0 7 4- f (Y Z ) 


(D.IO) 
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^ 3 T" 


0:X.rs 


^ J h + trr ^ "1 


El 


(B.ll) 


iM. ih : filial condition 

4(0) - 


.0 ^ 2.^ 


ii; derived from 


ff 

a, (&(0). * ^o'''''-i<l„(-c) {S-12) 


3 ju-j ^ • -0 

• wm> i^'* 


with initial oonclition 

(U) = 

J Nr r^) is obtained by 

,hc nrat oruor appro.iBaiion - C.. s) 

^ fT z ) as given m 

. iv,« corroD ponding terms to (X^. ,4.^^ 

adUinc. i<he con . -Pnr< a , P^ and ai 

„ ,6S The al>ova editions ior do- Eo 

equation (4.1 )• expansions of f and £, 

me trior series exp^xioj- 
aro derived by using Tay 

about ^oCOl 
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